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I. INTRODUCTION 


The Advanced Transport Operating Systems (ATOPS) program conducted by the 
Langley Research Center of the National Aeronautics and Space Administration is a re- 
search and development program aimed at developing capabilities for increased termi- 
nal area capacity, safe and accurate flight in adverse weather conditions including shear 
winds, the avoidance of wake vortices and reduced fuel consumption. Advances in mod- 
ern control design techniques and increased capabilities of digital flight computers coupled 
with accurate guidance information from the Microwave Landing System (MLS) make the 
achievement of some of these goals feasible. The development of the stochastic feedfor- 
ward/feedback control design methodology and its application to the design of a digital 
automatic landing system for a small transport jet aircraft was performed within the con- 
text of the ATOPS program. 

The main objective of a control system may be described as to enable the plant to 
track a desired trajectory, usually selected out of a given class of trajectories, as closely as 
possible in the presence of random and deterministic disturbances and despite uncertainties 
about the plant. Thus, a control system generally has a feedforward controller which tries 
to track a desired trajectory, and a feedback controller which tries to maintain the plant 
state near the desired trajectory in the presence of disturbances and system uncertainties. 
It is essential that the feedback law produce a closed-loop system which is stable about 
the desired trajectory so that small disturbances can be accommodated while maintaining 
the plant state near the desired trajectory once the feedforward control law has brought 
the plant state to the desired trajectory. 

The part of a control law which uses only the desired or commanded trajectory ex- 
plicitly will be referred to as the feedforward control law. On the other hand, the part of a 
control law which explicitly uses only measurements of the plant state will be referred to 


as the feedback control law. Even though in some cases involving nonlinear control laws 
the distinction between the feedforward and feedback control laws may become somewhat 
ambiguous, in linear control laws, the distinction is rather straight-forward. 

The design of the feedback controller has received considerable attention in the modern 
control literature in the last two decades; e.g., see [l] - [5], and the references therein. 
On the other hand, the design of the feedforward controller has received relatively little 
attention [6] - [10]. In this study, a combined stochastic feedforward and feedback control 
design methodology is developed and is applied to the design of a digital automatic landing 
system for the ATOPS Research Vehicle, a Boeing 737-100 aircraft. 

The feedforward control problem is formulated as a stochastic optimization problem 
and is imbedded into a stochastic output feedback problem [11], [12] where the plant 
contains unstable and uncontrollable modes. As the standard output feedback algorithm 
requires an initial gain which stabilizes the plant, a new algorithm is developed to obtain 
the feedforward control gains. The necessary conditions are shown to result in coupled 
linear matrix equations, implying that when a solution exists, it is indeed the globally 
optimal control gain. 

The formulation of the feedforward problem in a stochastic, rather than the standard 
deterministic, setting is significant in two ways. First, the class of desired trajectories 
from which the actually commanded path is selected can be effectively described as a 
random process generated by a dynamical system driven by a white noise process. The 
second, and more important, implication of a stochastic optimization formulation is the 
tacit understanding that “perfect tracking” is often not possible due to various reasons 
including uncertainties about, or vanadiums in, unc piaui paiamcucio, ww ^/icociu.c ui i^iaiio 
nonlinearities and unmatched initial conditions. Thus, questions about the robustness and 
sensitivity of the feedforward controller arise naturally in this context. 

A combined stochastic feedforward/feedback control methodology is developed where 
the main objectives of the feedforward and feedback control laws are clearly seen. Fur- 
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thermore, the inclusion of error integral feedback, dynamic compensation, rate command 
control structure, etc. is an integral element of the methodology. Another advantage 
of the methodology is the flexibility that a variety of feedback control design techniques 
with arbitrary structures may be employed to obtain the feedback controller; these in- 
clude stochastic output feedback, multi-configuration control, decentralized control [13] or 
frequency and classical control methods. 

Finally, a specific incremental implementation is recommended for the combined feed- 
forward/feedback controller. Some advantages of this digital implementation are the sim- 
plicity of implementation, the fact that trim values are not needed and that problems 
such as integrator wind-up can be largely avoided. The closed-loop eigenvalues using this 
implementation are shown to contain the designed closed-loop eigenvalues which would 
result if an incremental implementation were not used. 

A digital automatic landing system for the ATOPS Research Vehicle (a Boeing 737- 
100) is designed using the stochastic feedforward controller and stochastic output feedback. 
The system control modes include localizer and glideslope capture, localizer and glideslope 
track, crab, decrab and flare. Using the recommended incremental implementation, the 
control laws are simulated on a digital computer and interfaced with a nonlinear digital 
simulation of the aircraft and its systems. 
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n. STOCHASTIC FEEDFORWARD/OUTPUT FEEDBACK 

CONTROL DESIGN 


A control system generally contains a feedforward and a feedback control subsystem. 
The feedforward controller tries to track a desired input (or commanded) trajectory, while 
the feedback controller tries to maintain the plant state near the desired trajectory despite 
the presence of disturbances, random noises and system uncertainties usually by using 
error feedback. In the modern control literature, the design of the feedback controller 
has received considerable attention, while the design of the feedforward controller and its 
relationship to the feedback controls has received relatively little attention. On the other 
hand, classical control techniques have treated the design of both feedback and feedforward 
controllers jointly. 

In this section, we will formulate the design of a feedforward controller as a stochastic 
optimization problem. We will present the solution to this problem for two important 
special cases. Then we will present a control design methodology which combines the 
feedforward and feedback control designs and addresses various questions which arise in 
practical control law design. 

A. Feedforward Control - A stochastic Formulation 

In general terms, the objective of a control law is to enable the plant to track a 
“desired or commanded trajectory” as closely as possible, in the presence of disturbances 
and despite uncertainties about the plant. 

In this study, for purposes of discussing terminology, the part of a control law which 
uses only the desired trajectory, or the command state, explicitly will be referred to as the 
feedforward control law. On the other hand, the portion of the control law which explicitly 
uses measurements of the plant state will be referred to as the feedback control law. Even 
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though in some cases it is difficult to separate the feedforward and feedback controllers, in 
most linear control laws, the distinction is relatively straight-forward. 

Command Model. 

Consider the desired trajectory or command model 


Zk+l = <f>zZk + (k, ( 1 ) 

where Zk is a n, -component command state vector; the order, n z , of the command model 
is arbitrary, and can be higher or lower than the order of the plant to be controlled. Note 
that not all the components of the command state, Zk , need correspond to actual physical 
quantities such as plant states or measured sensor outputs. Let H z be a n 0 x n z matrix. We 
assume that only the vector H x Zk will be directly commanded as the desired trajectory. 
For example, in an altitude control law, only the altitude may be commanded as the desired 
trajectory; however, the command state vector may be defined with two components (i.e., 
n z = 2) possibly corresponding to commanded altitude and sink rate. 

From Eq. (1), it is clear that for an arbitrary command state history { Zk , k > 0}, it 
is always possible to determine the “forcing function” (or vector) {&, k > 0} such that 
the command model Eq. (l) holds. This can simply be achieved by solving for ft. If 
the command state history {«*, k > 0} is fixed or completely known a ‘priori’, then a 
control sequence which makes the plant track this trajectory can be obtained provided 
that the trajectory is realizable [5]. However, this control sequence would be a fixed, open 
loop control sequence corresponding only to that specific trajectory, rather than being a 
feedforward control law. 

In most applications, we are interested in designing a feedforward control law which 
can track any one trajectory selected from a given class of command trajectories, say z{£,). 
One method of specifying a class of trajectories is to specify a dynamical model driven by 
a random process with given statistics. 
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For example, consider the command model in Eq. (l), where the sequence {ft, k > 0} 
is specified as a vector random sequence with white noise statistics. The class of command 
trajectories thus specified would be the family of command state histories {z k , k > 0} 
which can be generated by any realization of the random sequence {ft, k > 0}. On the 
other hand, {ft, k > 0} may be specified as a random sequence with colored noise statistics. 
If ft itself can be obtained as the output of a discrete dynamical system driven by white 
noise, then the command state may be augmented to obtain a new system of higher order 
driven by white noise but still having the form of Eq. (1). 

It is clear that a large family of command trajectory classes can be expressed by the 
model given in Eq. (l) by appropriate selection of the system order n M , transition matrix 
<f> z , and the covariance of {ft, k > 0}. In the remainder of this study, we will assume that 
the random sequence {ft, A: > 0} in Eq. (l) has white noise statistics, unless specified 
otherwise. 

Feedforward Control Optimization. 

Consider a linear plant model of the form 


a*+i = <t>x *k + r u «* + r, z k + r f ft + w k , (2) 

where x k is the n-component plant state vector, u* the r-component control vector, w k the 
plant noise process, <f> x the plant state transition matrix and T u the control effectiveness 
matrix. The vectors z k and ft are the command state vector and the command forcing 
vector, respectively. 

It should be noted that when the matrices r z and T f are selected to be zero, the plant 
state x k does not depend on the command trajectory, which is the usual case. However, 
it is often desirable to include in the plant model, states which describe the error, or the 
deviation from the command value. To accommodate design models of this type, it is of 
interest to include the command trajectory terms at this point in the formulation. 
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Let H x be a n a x n matrix. The objective is to obtain a feedforward control law so that 
the plant variables H x x k continuously track the commanded variables H z z* as closely as 
possible when {z*, k > 0} belongs to a given class of command trajectories, L. 

By its nature, a feedforward control law is intended not to modify the stability, noise 
attenuation and robustness properties that are already present in the plant model. These 
properties are generally obtained by appropriate design of the feedback control law. For 
the purpose of designing the feedforward control law, we will assume that the plant model 
already incorporates the feedback control law. Thus, in this section, Eq. (2) represents 
the closed-loop plant model where the feedback control law has satisfactorily achieved the 
desired closed-loop objectives. In particular, we will assume that the closed-loop system 
is stable; i.e., all the eigenvalues of the state transition matrix <f> x lie inside the unit circle. 

We will consider feedforward control laws of the form 


u k — Kz z k ~ ?k» (3) 

where K z and K s axe control gains to be selected in order to track the command tra- 
jectory as closely as possible. It is important to note that Eq. (3) is one of the simpler 
control structures that can be selected. More complex feedforward structures should be 
investigated to extend the results obtained in this study. 

In order to obtain a set of feedforward gains K z , and K s , it is desirable to select a 
criterion or objective function which describes the goals to be achieved, and then optimize 
this criterion. Since our goal is to track H z Zk , an obvious selection would be 

1 N 

J = «!So 2(A' + l) E £ (g - " H ' ** " H ’ **> (1) 

While more general yet quadratic objective functions can be selected, some properties 
of this expression may be noted. Since only the tracking error is penalized, if a control can 
achieve perfect tracking then it will optimize the criterion. Thus, the optimal control will 
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result in perfect tracking when that is possible with the form of control selected in Eq. (3). 
Otherwise, the optimal solution will minimize a quadratic function of the tracking error. 
It should be noted that, in this context, perfect tracking implies that almost all command 
trajectories in the class t can be tracked by the plant using the feedforward control law 
in Eq. (3). 

While the lack of control weighting allows perfect tracking, when possible, it may 
also result in more control activity than desired. To accommodate such cases, it is always 
possible to include a non-zero control weighting term. However, this inevitably results 
in less than perfect tracking; in such cases, it may be a better policy to change (e.g., to 
smooth) the commanded trajectory so that tracking it does not require as high a level of 
control activity. 

In the following, we will use the objective or cost function 

&)(:) w 

which does not explicitly contain a control weighting term. The case given in Eq. (4) is 
obtained when 


Qn = Hi H x> Q 12 = -Hi H x , (6a) 

Q21 — —Hi H x , Q 2 2 = HlH x . (6b) 


Tile stochastic feedforward optimization problem can now be posed as; Given the 
plant model in Eq. (2) and the command model in Eq. (1), find a feedforward control law 
of the form of Eq. (3) which minimizes the cost function J in Eq. (5) . 

When the plant model is augmented by the command model, the feedforward control 
optimization problem posed is seen to be a stochastic output feedback optimization prob- 
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lem [6] - [8], where the “feedback” vector contains only the co mm a n d variables z* and ft.* 
Thus, it would appear that the stochastic feedforward problem can be embedded in the 
stochastic output feedback problem. In fact, when the command model transition matrix, 
cfrzi is stable, the output feedback algorithm can be directly used to determine the optimal 
feedforward gains K z , K i . 

However, most realistic command models require the use of unstable, in particular 
marginally (un)stable, systems. For example, the command model for the usual case of a 
constant command would have an eigenvalue of unity. 

Zk + 1 = Zk + 0 ft, z c = const. (7) 

Therefore, for the stochastic feedforward problem, the assumption that the command 
system is stable is not a realistic hypothesis, and is of limited use. Since the command 
model is not controllable, it is not stabilizable except when it is already stable. Thus, 
it is not realistic to assume that the augmented system is output stabilizable. Since the 
output feedback algorithm [9], [12], requires a stabilizing gain, it cannot be used to obtain 
the optimal feedforward gain. The fact that, in most cases of interest, the command 
model is not output stabilizable produces the major difficulty in determining the optimal 
feedforward gains. 

In overcoming this problem to obtain an optimal feedforward control law, we will 
consider two cases. However, first we restate the necessary conditions for the optimal 
stochastic infinite-time output feedback problem [4], [9]. 

P(K ) = <f> T (K) P(K) <}>{K) + C t K t RKC + Q , (8) 

S(K) = <j>{K) S<f> T {K) +TKV K t T t + W , (9) 

*Note that ft can also be included in the augmented state and, hence, in the “feedback” 
vectors, as is shown in Case 2 which is discussed later. 
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P{K) K S(K) = r T P(K) <f> S(K) C T 


( 10 ) 


4>{K) = *-TKC , (11) 

P(K) = T t P(K)T + R , S(K) = CS(K)C T +V , (12) 

where P(K) and S{K) represent the discrete cost and state covariance matrices, respec- 
tively, when the gain K is used. Q and R are the discrete state and control weighting 

A A 

matrices, respectively; W and V are the plant and measurement noise covariance matrices, 
respectively; C is the output gain matrix, as described in [12], and K is in the stability set 
S. 

Case l. tt* = — K g zk- 

The augmented system for this problem can be expressed as 


( z *+ iN ) _ ( r,\ / x k \ 

\2k+l ) \ 0 <t>z) \Zk ) 

+ (o)“‘ + (“‘ + ft r!C ‘) 

(13) 

Vk = C M z k + v k = (0 

c -)(2) + ”* ■ 

(14) 

ttfc = -K a y k = 

— K z C z 2^ — K z Vk 

(15) 


In order to determine necessary conditions for the case considered, we rewrite the 
general necessary conditions given in Eqs. (8) - (12) after partitioning the matrix equations 
according to the partition in Eq. (13). This results in the following equations. 

P** = 4>lP*.4>* + Q.m , ( 16 ) 
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Pxz = 4>Z P x , <f>, + 4>Z Px*(r, - r x k z c z ) + q xz 




(17) 


S xz = <f> x s xz <f>Z + (r, - r x k z c z )s zz <f> T z + w xz , (is) 

&zz = <t>z &zz 4>Z H" W zz , (19) 

K z = P xx TZ[P xx <f>x S xz + (P xx T z + P xz <t> z ) S zz ] Cj S~ z , (20) 

where 


Pxz = P = Tl Pxx Tx + R 

s zz = s = c z s zz cZ + V 


( 21 ) 


( 22 ) 


From these equations, it may be noted that the optimal feedforward control gain 
matrix, K z , does not directly depend on P zz nor S xx . Furthermore, recall that, strictly 
speaking, the necessary conditions hold when the closed-loop system is output stabilizable. 
For the current problem this requires that <f> z be stable. 

It should be noted that when <j> z is unstable, the covariance of the command vector 
grows without bounds; so that Eq. (19) does not have a non-negative definite solution, 
S zz . On the other hand, Eq. (17) has a finite solution, P xz , under relatively unrestricted 
conditions [13]. In particular, if 


p{4>z)p(<l>Z)< 1, (23) 

then a finite P xz satisfies Eq. (17), where p(<f>) denotes the spectral radius of the matrix 
<f>. Since <f> x is the closed-loop plant transition matrix, and is assumed to be stable p[4>Z) 
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is strictly less than unity; if necessary, it could be designed to be smaller. On the other 
hand, in most cases of interest, a p{<t> z ) value of unity is sufficient to model the desired 
trajectory class. 

Similar comments apply to the solution of Eq. (18) since the matrices involved have the 
same eigenvalues; i.e., p{<f> x ) p{<f> z ) is equal to p{<j> z ) p{4> z ) . However, the forcing function 
in Eq. (18) contains S zz . When S zz grows without bounds, so will S xz . However, this 
does not necessarily imply that the optimal feedforward gain, K zt will also grow without 
bounds. 

To investigate the nature of the optimal solution in the limiting case where the out- 
put stabilizability condition does not hold, we rewrite Eq. (18) in the following manner. 
Assuming that S zz is nonsingular, let 


Szz — S zz S zz ; (24) 

then S xa satisfies the equation 

Sx, = [*. $„ + (r, - T x K z C a )][S„ <f>? S- 1 ] + W xz S- 1 (25) 

Further assuming that <f> z is nonsingular, and manipulating Eq. (20) results in 

s„ t T . s~' = - w„ S~‘] . (26) 

Substituting Eq. (26) into Eq. (25), and manipulating, we find that 


c rr_w e-ii-i 

r r ZZ “ZZ J 


A 


'f'Z 


-i C I (P 
— 'TX *XZ 


r* zr /-* \ 
x x *^z ^z) 


t ir To t ? r i — l i 

”xz[^zz — <Pz (**) 


Now, in many cases of interest, when <f> z approaches an unstable matrix, the covariance 
of the command state, S ZZy grows without bounds, while S zz l vanishes. In this report, we 
will limit attention to cases where the inverse of the covariance of the command state 
vanishes, unless stated otherwise. This results in simplifications in Eq. (27) and Eq. (20). 
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We will further assume that all the command state components are known without any 
error. This is certainly a reasonable assumption, since we cannot command a trajectory 
that we do not know; thus, we set the measurement noise covariance, V, to zero; i.e., 


V' = 0 , and C Z = J . (28) 

The necessary conditions for optimality for this case can now be expressed as: 

Pxx = ^Z P XX $x + Qxx > (29) 

Pxm = <t>Z Pxx + <t>Z Pxx{Tx - T x K z ) + Q xz , (30) 

S xz <f> z <f> x S xz = T z T x K z , (31) 

K z = P xx T-Z[P*X $2 "k Pxx r* + Pxx < t > x\ > (32) 

where P xx is given by Eq. (21). 


It should be noted that P xx is independent of the feedforward gain matrix, K z , but 
depends only on the (closed-loop) plant transition matrix, <j> x> and the weighting matrix, 
Q xx . Similarly, the covariance matrix, S zx , is also independent of the feedforward gain 
matrix, K z . However, S xz has been eliminated from the necessary conditions and does not 
explicitly appear in these equations. 

It is of interest to note that the solutions of Eq. (30) and Eq. (31), P xz and S xz , re- 
spectively, are linear functions of the feedforward gain, K z , as these are standard Lyapunov 
equations. Since P xx , P xx and 4> x are independent of K z , Eq. (32) is also linear in the 
gain K z . Thus, obtaining the optimal feedforward gain, K z , does not require the solution 
of highly nonlinear matrix equations, but can be obtained by solving a set of coupled but 
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linear equations. An algorithm to solve for the optimal feedforward control gains, K z , will 
be shown in a later section. 

CASE 2. u* = -K s (Sk + v k ). 

In the preceding case, the feedforward control law was restricted to using only the 
command vector, z k . In this section, we will consider the additional use of the command 
model forcing function (or vector), ft, in the feedforward law. Since ft is necessary to 
determine the succeeding command vector, z k +u it contains lead information and can 
play an important part in satisfactorily tracking the desired trajectory. 

The augmented system for this problem 


/ Xk+i ( 4>x r, — T s K m T s \ fx k \ /r*\ ( Wxk \ 

Z k +1 = ° <f> z I ) I + 0 tft + w,* » ( 33 ) 

V Ck+1 J V 0 0 0 J \Ck J \o J \w ik J 


yk = (k + vk = {oo i) 



+ Vk 


(34) 


Partitioning the necessary conditions given by Eq. (8) - Eq. (12) according to the 


partition of Eq. (33), results in the necessary conditions for the problem considered here. 


Pxx — 4 > 's Pxx d" Qxx (35) 

Pxz = 4>l Pxx <t>z + <t?z Pxx (r, - r x Kz) + Qzz (36) 

= Pxx r l\Pxx r f + Pxz] W s [W { + V;}- 1 (37) 

Pxx = P xx T x + R (38) 
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It is important to note that the solutions P xx and P xz to Eq. (35) and Eq. (36), 
respectively, are precisely the same as the solutions to Eq. (29) and Eq. (30) in Case 1. 
Thus, if P xx and P xz are computed when obtaining the optimal gain for the command 
vector feedforward, the same matrices can be used in obtaining the optimal gain for the 
forcing vector feedforward. 

Furthermore, it should be noted that the command forcing vector feedforward gain, 
K i , given by Eq. (35) - Eq. (38) is optimized for an arbitrary command vector feedforward 
gain, K z . In other words, Eq. (35) - Eq. (38) hold for an arbitrary gain, K z , not only for 
the optimal K z . 

Finally, note the simplicity of the expression for K s . Once the weighting matrices have 
been appropriately selected, and the feedback gains have been obtained so that satisfactory 
feedback characteristics are achieved, it is possible to compute K i on-line using the current 
values of the closed-loop plant transition matrix, <f> x , and the control effectiveness matrix, 
r*. Thus, the simplicity of the necessary conditions for the optimal make it usable as a 
feedforward gain-scheduled controller. Similar comments apply to the optimal feedforward 
gain, K z , for the command vector. However, in the latter case, the necessary conditions, 
although linear, are not as simple and easily computed as for K$. 

It is of interest to consider the case where the command forcing vector is known with 
no error (i.e., V s =0), so that & can be fed forward without noise. In this case, Eq. (37) 
becomes 


K S = P xx T x [Px X Ts + P XZ } • (39) 

Thus, the optimal when ft is perfectly known at the k th sampling instant is 
independent of the covariance, W$, of ft. From an alternate point of view, the gain K s can 
be decreased, or altered, by appropriate selection of the “measurement noise” associated 
with the forcing vector. 
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An Algorithm for K z . 

As mentioned earlier, the necessary conditions for Case 1 which considers the feed- 
forward of the command vector, z*, are linear functions of the gain K z , which can be 
seen by observation of Eq. (29) - Eq. (32). Although linear, these matrix equations are 
coupled; so that an explicit expression for K z cannot be easily obtained. In contrast, the 
gain K t , the feedforward of the command forcing vector, can be easily solved for in an 
explicit expression as shown in Eq. (35) - Eq. (38) . 

In order to develop an algorithm which results in the optimal feedforward gain, K z , we 
will make use of a basic principle of linear operators. Let l be a linear transformation from 
some p-dimensional real linear space into itself, and let {e,-, 1 < * < p} be an arbitrary 
set of basis vectors spanning the space. To an arbitrary vector K in the space, having the 
representation 


K=Y1 K ' e < - Ki* R 

associate the column vector KeR p 


p 

E 

t=l 


(40) 


K = 


r*i\ 

K 2 

\Kj 


(41) 


Then, the matrix representation, L, of the linear transformation, l, satisfies the equation 


£(ey) — Lj — ^ Li] ti , 1 < ] < P 

»=i 


L= (Li|L 2 |...|L p ) 


(42) 


(43) 


If K is given by Eq. (40), then 


e(K) =LK 


(44) 
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Now, rewrite Eq. (32) in the form 


r^-P XX ^X 5 Z * + Pxx r*x 4* Pxx <f>z ] P XX -K* — 0 


(45) 


Equation (45) can be viewed as a linear transformation from the space of rx n z matrices 
into itself. We will view Eq. (45) as corresponding to 


1{K Z ) +L 0 = 0 

where L 0 is a r x n z matrix. 

Consider the basis {ey, 1 < j < rn z } defined by 


Cl = 


/I 0 .. 

0 0 .. 

• • • 

• • * 
• • 

Vo 0 .. 

Cr+1 = 


Of course, other basis selections are also possible. 
From the preceding discussion, it follows that 


0\ 



o 

O 

o 

0 



1 0 ... 0 

I 

J 

«2 = 

• • • • 

• • « t 

• * • • 

oj 



Vo 0 ... 0 J 

f° 

1 ... 0\ 



0 

0 ... o' 



: 

I * • ; 

» 

etc. 

Vo 

0 ... 0 ) 




(46) 


(47) 


(48) 


LK g + L o = 0, 


(49) 


K z = —IT 1 L 0 , 

whenever L is invertible. This leads to the following algorithm. 

Algorithm: 

STEP 1. Solve the Lyapunov equations 


(50) 


Pxzj — 4>Z Pzzj 4>M P ** e i » 1 ^ j < rrl z 


(51) 
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S X zj 4*z <f>xSxzj — r* Cj , 


1 < j < rn z 


(52) 


STEP 2. Compute L 


Lj — xx 4*x Sxzj "H P tzj 4*z\ P xx Cy (53) 

STEP 3. Form the L 

L=[L 1 \L 2 \...\L rnM ). (54) 

STEP 4. Solve the Lyapunov equations 

P XZO = 4*x ^ 1*0 "P 4x PxX P* Qxz (55) 


Sxzo 4>z 4x Sxzo — Pi 


(56) 


STEP 5. Compute L 0 and form L 0 


Ijq — Tx\Pzx <t>x Sxzo "t" P xx r z + Pxzo 4* 


(57) 


STEP 6. Solve the linear equation 


LK g + L 0 = 0 


(58) 


orpTHTJ fr p n .. . 1/ J aL- ( J! a) - a..a 

w X *1.* 7 • a orm 71 jg unu wiiiyuic tliC gi dUlCilt OS 41 bCSl 


Pxz = K Pxx 4>z + 4>l Pxxpz - T x Kx) + Qxz 


(59) 


Sxz 4>z ~ <t>x Sxz — r z r z k z 


(60) 


18 



( 61 ) 


Px*k m r^" \p x x dx $xz + Pxx r* + P X z 4 > z\ || 

STEP 8. Compute the forcing vector gain, K s . 

K s = p- 1 r l\p xx r f + P xz ] w s [W s + v s f 1 (62) 

The algorithm uses the fact that S xz and P xz , although ‘linear’, are not homogeneous 
functions of K z \ so that the non-homogeneous part is separated in the algorithm and 
combined in a single term, namely L 0 . 

It should be noted that the algorithm is not iterative, so that convergence questions 
do not arise. The problem solved is a set of coupled linear equations, and the solution is 
functionally obtained by inverting a rn z x rn z matrix, namely L. 

On the other hand, the algorithm requires the solution of 2(rn z -f 1) Lyapunov equa- 
tions and the inversion of a square matrix of dimension rn z . When dealing with ill- 
conditioned matrices and high order problems, the accuracy of the solution may be of 
concern. Usual techniques to improve accuracy may be used in such, as yet hypothetical 
cases. Importantly, it is possible to test the accuracy of the solution, K zt by computing 
the ‘gradient’ of the cost function with the exception that S zz has been eliminated from 
the expression for the gradient. The norm of this pseudo-gradient which is computed in 
Step 7 of the algorithm provides an indication of the accuracy of the solution obtained. 

The solution is unique if L is invertible. In this case, the feedforward gain obtained 
is the global solution to the stochastic optimization problem considered. It is desirable to 
obtain conditions which determine the invertibility of the matrix L in terms of parameters 
such as <j> x , 4> z , r x , etc. which are more directly related to the physical aspects of the 
control problem under consideration. On the other hand, the author’s experience on this 
problem which has been necessarily limited, has always resulted in an invertible L, hence 
a unique solution of the global optimal solution. 
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Finally, it may be of interest to note an identity which can be used in the algorithm 
as an alternate expression for Eq. (53), Eq. (57) and Eq. (59). 

rx 4*x $xz Prz r, + P X z 4>z\ Pez K z 

— {p*z[^z $xz + r* — r z k z \ + Pxz — R k z (63) 

= r l (P« S XM + Pxz) <f>z — R Kg . (64) 

where Eq. (31) has been used in Eq. (63). The expressions in Eq. (63) and Eq. (64) may 
be interchanged when desirable. The latter expression is somewhat simpler, particularly 
when the control weighting matrix, R , vanishes. 

Constant Command, <f> z = I. 

The most straightforward way in which a constant command can be modeled would 
be to simply select a command model where 

1 = *k + ft • (65) 

As a constant command is a commonly used desired trajectory, we will investigate this 
special case in more detail to obtain the feedforward control law for it. As will be seen, 
it is possible to obtain an explicit expression for the feedforward gain, K z , for this special 
case. 

The necessary conditions in Eq. (30) and Eq. (31) can now be solved explicitly for 

x\ z bU ujbaiu 


Pxz = (I- PxxlTz - Tx Kz) + Qxz] 


( 66 ) 


s** = (/-<M“ 1 (r*-r I .K',) 


(67) 
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Substituting these expressions into the necessary conditions expressed in Eq. (64) and 
manipulating, we obtain 


K a = (T? p xx T x + R)~ l Tl[ Pxx T x + (/ - tf)- 1 Q xz \ , (68) 

Pxx = Pxx(I - tx)"' + (/ - <£)- '</>* Pxx , ( 69 ) 

Pxx = Pxx + Pxx M* - <t>x)~ X + (/ - +l)- X *Z Pxx , ( 70 ) 

where P xx is given by Eq. (29). 

It should be noted that since the closed-loop plant transition matrix, <f> x is stable, all 
its poles are strictly inside the unit circle. Therefore, ( I — <f> x ) is invertible. In this case, 
the existence and uniqueness of an optimal solution is determined by the singularity of the 
matrix (Tj p xx T x + R). 

For the constant command case considered here, it is therefore not necessary to use 
the algorithm given in the previous section. Finally, it is possible to obtain a similar 
expression for the command forcing vector feedforward gain, if f , by simply substituting 
Eq. (66) into Eq. (37). 

A Simple Example. 

To illustrate the stochastic feedforward control law optimization developed in the 
preceding sections, we will consider a simple 1** order example with a constant command 
model. 


Xk+i = <P x k + 7 u k + w k , \<p\ < 1 , (71) 

where is a white noise sequence with zero mean uncorrelated to the command state 

and forcing vector, z k and &, respectively. We model the class of desired trajectories by 
the l 4t order command model 
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Zk+ 1 = Zk+Ck 


( 72 ) 


where we model {&} as a zero mean white noise sequence uncorrelated to the initial 
command z Q which may have a non-zero mean. 

Suppose that we would like h x Xk to track h z Zk as closely as possible at all sampling 
instants irrespective of the control activity required. Now, using the results of the last 
section, we have 


Q»g(l ~H < p) Qzz hg 

(1 — SP 2 )(1 — Sp) (1 - <p) 2 ~ (1 - <p) 2 


2 Qzz 2 9n 2i 

_ (l-y?) (1 - y?) 2 _ (1 ~ <p) Qx Z 2l 

* 1 2 Qxx mT} 7 Qxx 1 

where R has been assumed zero. If we further substitute —h x h z for Q XM , 
null -7,, 


( 73 ) 

( 74 ) 

h x h z for Q xx , and 


K x = 


(1 - ip) h z 


1 h s 

Using the optimal feedforward of the command vector, Zk, 


( 75 ) 


Xk+l =<px k + (l-<p)-—Zk + Wk 

n x 

To see the effect of this control law, suppose that at some sampling instant k, 


( 76 ) 


hx x k — h z z k 


( 77 ) 


then 


h z tt \ h x h x 

Xk + 1 = <p-£- Zk + (1 - lp)^- Zk + Wk = Jj- Zk + Wk 


( 78 ) 
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It is seen that, as would be expected from an feedforward controller designed for 
a constant command, this control law maintains x k+i at a constant desired value (with 
the exception of plant noise effects which are attenuated by the feedback design). This 
feedforward controller drives the state to the desired trajectory, based on the assumption 
that the command vector is most likely to remain constant. Note that 


z k+ i = E(z k+ i\z k ) = z k , (79) 

so that the least-squares estimate of the next command h z z k+i is in fact the last command 
h z z k . 

Thus, when the feedforward controller is limited to using only the command vector, 
based on a constant command model, the design performs exactly as would be expected, 
and drives the state to the desired constant, and then maintains the commanded value. 
If the command state, z k , is not completely constant, but moves slowly, the state will 
track the movement with some error. It may be noted that the Command-Generator- 
Tracker (CGT) deterministic feedforward controller for this problem is the same as the 
one obtained here; i.e., Eq. (75). It should also be noted that, for this example, non-zero 
control weighting (i.e., R > 0) produces less than perfect tracking. 

Now we remove the restriction that only the present command value, z k , be used in 
the feedforward law; by allowing the current value of ft to be used. Assuming that ft is 
perfectly known at sample k, we set V" f equal to zero. Using Eq. (39), the optimal gain for 
the forcing vector, ft, can be found to be 


K' = - 


<P 


K z + 


(1 + V?) h 2 


[(I-*?) 7 h x 

Substituting the optimal K z from Eq. (75) into Eq. (80), and manipulating 


(80) 


K i=~ 


1 h* 
7 K 


(81) 
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The full feedforward control law now becomes 


(l <p ) h z 1 h z 

Uk = - T- Zk + - 7— fjfe 

'I h x 7 h x 


( 82 ) 


/i h 

**+1 = Zfc + (1 - <p) Zk + ft + tu* 

/lx /lx 


(83) 


Assume, as before, that at sample k perfect tracking was achieved; i.e., Eq. (77) holds. 
Then, 


zjfc+i = (** + ft) +Wk = J t 1 Zk+1 + Wk • ( 84 ) 

When the feedforward law is not restricted to Zk alone, it is seen that the state 
tracks the command perfectly, neglecting the plant noise effects. This improvement in the 
tracking performance is due to the availability of accurate lead information in f*. It is 
clear that, in this simple example, the feedforward law given by Eq. (82) will track the 
desired trajectory at every sampling instant, as long as the plant parameters <p and 7 are 
perfectly known and do not vary. 

Although illustrative of many true trends, this example is of low order and does not 
represent all the complexities of a realistic design problem. 

B. A Stochastic Feedforward/ Output Feedback Design Methodology 

Having both feedforward and feedback control design techniques available for use, it 
is necessary to also have a methodology which combines these two control designs to solve 
realistic and practical control design problems. In particular, the control design method- 
ology should have the capability to: 1) accommodate a variety of control law structures, 
2) allow the use of different control design techniques to achieve desirable characteristics, 
and 3) have a simple and practical implementation free of the often-encountered problems. 

In many control design problems, the specification of the control objectives also imply 
and sometimes require a particular structure for the control law. For example, if it is 


24 



necessary to have zero steady state tracking error in some variable in response to a constant 
command despite small variations in the plant parameters, then an error integral feedback 
structure is necessary to achieve this objective. Whereas if a bias error in that variable 
can be tolerated, then only an error feedback is sufficient.* 

Similarly, in many problems the unavailability of sensors which accurately measure 
every state variable makes it desirable to use dynamic compensation in the feedback loop. 
In other cases, a decentralized control structure may be desirable. A methodology which 
can accommodate a rich collection of the combined feedforward and feedback control law 
structures is desirable. 

On the other hand, a large variety of feedback control design techniques are now avail- 
able. These include modern control design methods such as stochastic output feedback 
or full state feedback techniques, multi-configuration control (MCC) techniques, decen- 
tralized control techniques [13], as well as classical control and frequency domain design 
techniques. As most of these design methods focus on achieving certain desirable control 
characteristics, it is desirable to have a methodology which accommodates the use of many 
design methods. 

Finally, many modern control design techniques sometimes fail to adequately consider 
the digital implementation of the control law. For example, methods of avoiding the use 
of trim values, eliminating integrator wind up, minimizing the effects of control rate and 
control position limits, asynchronous sampling of sensor outputs, delays, etc. ought to 
be integral concerns of the design methodology. A design approach in which as many as 
possible of the practical implementation issues are explicitly treated is desirable. 

To develop a methodology which formulates a combined feedforward and feedback 
control design problem which addresses most of the realistic issues raised above, consider 
a linearized plant model of the form 

Sfc+l = <l>x 2k + T* Uk + w x k + d* , (85) 
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where the state Xk, the control u*, correspond to the total values of the corresponding 
quantities rather than being perturbations about their trim values, tu* is the white plant 
noise and d x is a constant vector arising from the linearization of the nonlinear plant about 
some operating point. When the plant operating point changes, all the plant parameters 
such as <f > x , r*, r*, r f and d x may change; however, it is assumed here that the changes in 
these parameters occur at a much slower rate than the variations in the state Xk , controls 
ujt, etc., so that the plant parameters are assumed constant except when explicitly stated 
otherwise. 

The feedback vector is assumed to admit a linearization of the form 

Vk = C x £ k + v k + d y , (86) 

where d y again is a constant vector which depends on the point about which the lineariza- 
tion is obtained, and Vk represents the measurement noise assumed to be a white noise 
sequence. 

As in the previous section, we will limit attention to the class of command models 
which admit a discrete stochastic model of the form 

*fc+i = <t>z Zk + ft > (87) 

where {ft , k > 0} is a sequence of uncorrelated random vectors with zero mean, such that 


E(tk = 0 

£(ft Vk) = 0 

(88) 

£(ft<)=0 , 

o 

II 

o 

(89) 


We consider the basic control objective to be the design of a combined feedfor- 
ward/feedback control law where only the vector, y*, is used in the feedback law, and 
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H x %k tracks H z Zk as closely as possible for a given class of commanded trajectories in the 
presence of disturbances and despite uncertainties about the plant. 

In many applications, the designer may be interested in a control law structure which 
contains dynamic compensation and integral feedback to achieve certain objectives. Fur- 
thermore, it may be desirable to have an inner loop/outer loop structure or to have some 
fixed (previously designed) filter provide some state estimates. Such dynamic subsystems 
as an inner loop or a filter can be included in the plant model described in Eq. (85). 

Furthermore, in many cases it is desirable and important to distinguish between the 
variable describing the commanded control value and the actual position of the control 
(the physical quantity such as an elevator position, or an electrical current). These two 
quantities, the commanded and the actual control values, are usually related through an 
actuator subsystem which should be modeled and included in the plant dynamics model 
in Eq. (85). 

In other cases, the designer may decide on a control rate command structure rather 
than using control position commands. This has a variety of advantages such as low-pass 
filtering the control command thus reducing unnecessary control activity, and providing the 
commanded control position for feedback when the actual control position is not measured 
and used for feedback. 

In most of these cases, it is necessary to use a formulation which accommodates 
output feedback for the controller. Here, it is assumed that the types of control structures 
mentioned above have been included in the plant dynamics and the feedback vector models 
given by Eq. (85) and Eq. (86), respectively. 

Error Formulation. 

While the plant and feedback vector models can be used in the form given by Eq. (85) 
and Eq. (86), it is often convenient and desirable to have state variables which represent 
the error in some of these variables instead of their values with respect to some fixed frame 
of reference. In cases where a plant state variable is being commanded by a command 
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model state variable, defining the error as the difference is straight forward. Also state 
variables corresponding to the time rate of change of already defined error terms can be 
easily obtained. However, since many variables are not directly commanded, only some of 
the state variable errors terms can readily be defined in this manner, while the remaining 
ones maintain their previous representations. 

Let H be a n x n z dimensional matrix relating the command model state, z*, to the 
plant model state, 2*, and let 


Xk = x k - Hz k (90) 

Substituting Eq. (85) and Eq. (87) into Eq. (90) and manipulating, the dynamics 
equations for the new state vector, Xk, can be obtained. 

Xk+i = <j> x Xk + TiUfc + (<f> x H — H </> 9 )zk — H $k + tu*fc + d z . (91) 

It is seen that while the plant model in Eq. (85), as is usually the case, does not 
depend on the commanded model state and forcing vectors, the description of the plant in 
terms of errors introduces such terms; hence, the inclusion of command state and forcing 
vector terms in the previous section dealing with the design of feedforward control. It 
should also be noted that the control vector u* in the plant representations given by Eq. 
(85) and Eq. (91) is the same. So that the new plant representation still deals with the 
same control values that drive the plant. 

Now consider the feedback vector in Eq. (86). The components of the feedback gen- 
erally correspond to sensor outputs, filter outputs, dynamic compensation states, control 
command values, etc. It is possible to use the vector y* by substituting Eq. (90). 

yk = C x (xk + H Zk) + Vk + d y = C x Xk + C x H Zk + Vk + d y (92) 

On the other hand, it is often convenient to think in terms of error feedback. So that 
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a new feedback vector consisting of the error terms defined in Eq. (90) may also be used 
in the form of 9k — C X H z k \ i.e., 

9k - C x Hz k = C x x k +v k + d v (93) 

Since the command state, z k , is known at the k th sampling instant, implementation of the 
error feedback vector in real time is clearly possible. 

To maintain further generality, we will consider the plant and feedback models in the 
form given below 


^*+1 — x k + r* u k + r M z k + r f $ k + w xk + d x (94) 

y k = C x x k + C M z k + Vk + d y , (95) 

where the standard error feedback case shown above would correspond to C K being null. 

From the limited experience and experimentation performed in this study, it appears 
that better performance is achieved when the error formulation is used in as many vari- 
ables as applicable. Due to the limited time available, the reasons for the differences in 
performance or methods for best selection of H were not investigated in this study. 

Dynamic Compensation and Integral Feedback. 

In many cases, it is of interest to include dynamic compensation in the control law to 
achieve particular objectives. The objective may be to estimate a variable for purposes of 
feedback, or to provide more robustness or insensitivity. 

Similarly, it is usually of interest to have integral feedback of the tracking error, or of 
equivalent variables, in order to obtain a type 1 system. To accommodate these often used 
control structures, augment the plant state model given in Eq. (94) by the compensator 
and integral error models 
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Cfc+1 = <f>c Ck + U c k + Wck , (96) 

Ik +1 = Ik + A t(H v y k - Hz Zk ) + tvik , (97) 

H x = H y C x . (98) 

Note that the “integrator” is a digitally implementable accumulator. Also note that 
the tracking variables, H y yk, have been assumed to be a linear combination of known or 
measured, yk, rather than possible unavailable state variables in Xk- 

The dynamic compensator is assumed to be of order n c which can be selected arbi- 
trarily, according to the desired objectives. The compensator state transition matrix, <f> c , is 
also arbitrary, and should be selected in accordance to the cost function which will produce 
the closed-loop compensator. The white noise sequences w c k and wik are included largely 
for generality. They could be interpreted as round-off error, variations in H y , jitter, etc.; 
however, most importantly, they can be used as design parameters which modulate the 
optimal feedback gains. 

The basic form of the combined feedforward/feedback control law is assumed to be 
u x k — — K xy yk H xc Ck K x i Ik H xz Zk K X f -I- u x , (99) 

u c k = —H C y yk — K cc Ck — K c i Ik — K cz Zk — K ci Ck + « c • (100) 

Note that, for the augmented problem, the feedback vector, y*, is also augmented by 
the compensator state, Ck, and the error integral variables, J*. The feedforward control 
law has been constrained to use only the command state and forcing vectors, Zk and £&, 
respectively. Finally, tt x and ti c are unknown constant (with respect to k) variables arising 
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from the fact that the variables used are not the perturbed values, but correspond to the 
total values including the trim values of the variables. As an incremental implementation 
will be used, it is not necessary to actually compute the constant vectors. 

It should also be noted that the control vector, may be selected to be the rate of the 
actual control position variables by augmenting the plant model accordingly. This would 
result in a control rate command structure. 

Feedback Design Model. 

Suppose that the feedforward control sequences {u* fc , k > 0} and {u * k , Jfc > 0} pro- 
duce a satisfactory trajectory. Then, when no plant noise or measurement noise is present, 
then trajectory will be given by 


*k+i =<f>xX* k + T x u* xk +T x z k +T i ft + d* , (101) 

c *+l = &C c k "l" u ck > (102) 

n + 1 = Ik + A t(H v y* k - H , Zk) , (103) 

y* = C x x* k + C z z k + dy . (104) 


Since the actual plant, compensator and integral states evolve according to Eq. (94), 
Eq. (96) and Eq. (97), respectively, the deviations, or error, in these variables can be 
defined as 


X k — Xfc — X k t U x k — t l x k u xk > U ck (105) 

yfc = y* - y* , h = h - » c k = c k -c* k . (106) 
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Manipulating, it is seen that the error in the state has the dynamics 


Xk+l = <f>k Xk + T x + w xk > 

(107) 

Ck + 1 = <t>c Ck + U c k + W c k , 

(108) 

h + i = h + At H x x k + (ttfjfc + At H y v*) , 

(109) 

V k — C x £k + t/fc 

(110) 


The deviations in the actual state relative to the desired trajectory are seen to be due 
to plant and measurement noise, and initial condition mismatch. Of course, in practice 
these deviations are also due to changing plant parameter values, unmodeled nonlinearities, 
unmodeled dynamics, sampling errors, etc. 

Since all the terms containing the command model have canceled, the deviation about 
the desired trajectory is seen to be independent of the command state. Thus, the possibly 
unstable command model has no direct impact on the feedback control law design. Where 
highly nonlinear effects which involve the command state exist, the command may not 
cancel; however, this is not a usually encountered case. 

Therefore, the design of the feedback control law can be done largely independently of 
the feedforward control law. The usual major objectives of feedback, such as stability about 
the desired trajectory in the presence of disturbances and despite uncertainties about, and 
variations in, the plant models, cam thus be pursued using the dynamical system describing 
the deviations about the desired trajectory in Eq. (107) - Eq. (110 ). 

The feedback vector, in this case, is taken to be (yj cj /J) . So that a control law 
of the form 
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«xfc = -K xy y k - K xe c k - K xI Ik 


( 111 ) 


ttcfc — KcyVk K ce C k K e i I k (H2) 

can be designed using any modem or classical control design technique. In particular, 
note that stochastic output feedback [11], [12], multi-configuration control or decentralized 
control techniques [13] can be used for this purpose. In the following, it will be assumed 
that the feedback control law thus designed stabilizes the closed-loop system. 

Feedforward Control Model. 

Having designed a satisfactory feedback control law, recall that the feedforward control 
in Eq. (101) - Eq. (104) is arbitrary, and can now be selected using a stable closed-loop 
system. 

Consider the following change of variables in the control vectors u* xk and u* fc . 


H' xk = < k + K xy y k + K xe cl + K xI II (113) 

= <k + Key y k + Kee cl + Kel I* k (114) 

Substituting these expressions into Eq. (101) - Eq. (104), we obtain the feedforward 
control model 


4+1 = (^* - r* K xy c x )xi - t x K xe cl - r x k x1 It + r x n\ k 

+ (r, - r x K ty c g )z k + r s $ k + d* x (115) 

cj + i = —Key C x x* k + (^c — Kce) c k — Kel I k + Me* — K C y C M Z k + d e (116) 
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(117) 


n+1 =Ik+M H x xl - (H z - At Hy C z ) z k + d } 


Zk+i =</>zZk + ik (118) 

where 


dl=d X ~ T X K xy dy , d* = -K C y dy , d} = At Hy dy , (ll9) 

are constant vectors depending on the trim conditions for the particular operating points 
used. 

From the discussion above, it is clear that the feedforward control design problem is 
one finding control sequences {jx* fc , k > 0} and , k > 0} in terms of a given subset of 
the augmented state variables z£, c£, /£, and ft, when the command state transition 
matrix <f> z , is not necessarily stable, such that {H y y£ — H z z* , k > 0} is as “small” as 
possible. 

A full analysis of the many interesting cases where different subsets of the augmented 
feedforward model state are selected is beyond the scope of the current study. Only the 
two cases solved in the preceding section will be treated in some detail. However, the full 
state feedback case is worthy of note. 

Consider the case where the feedforward control law form is unrestricted and {& , 
k > 0} is a Gaussian white noise process, with all initial conditions also being jointly 
Gaussian. With a quadratic cost function, it is well-known [4] that the optimal control 
is the solution to the LQG problem. The case where the plant contains unstable and 
uncontrollable modes, has been treated by the author as a disturbance accommodation 
problem [8], [7], It is clear that the most accurate feedforward control law would be 
obtained by this unconstrained solution. 
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While this solution has a variety of desirable characteristics, it also has a disadvantage. 
It requires the computation of the augmented desired trajectory x£, c£, J£, Zk, ft as a 
part of the feedforward controller. With the increasing speed and memory capabilities of 
flight computers, this is not necessarily impractical, but will not be pursued further in this 
study. 

Now consider the case where the feedforward control is restricted to the form 

n' xk = -K xz Zk - K xs ft - nl (120) 

Mcfc = ~Kcz Zk ~ K cs ft - n* (121) 

Using a quadratic cost function of the form of Eq. (5) which includes the most common 

case 

J nthrf £ l|/r « y '“ - H - - < 122 > 

it is possible to obtain the optimal feedforward gains K xx , K xi , K ez , and K c{ . 

It is important to note that the treatment of the feedforward control law in the 
preceding section accommodates the cases where any one of the feedforward gains is set 
to zero or some other constant. 

To obtain the total control law, recall that from Eq. (105), 


u xk — U x k + “*fc > u ck — U c k + * (123) 


ttxjfc — Vk K xc Kxi Ik d" (124) 


Uck = —Key Vk — Kcc Ck — K C I Ik + Mcfc (125) 
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The combined feedforward/feedback control law with dynamic compensation and in- 
tegral error feedback can be obtained by closing the loop on the compensation in Eq. (96) 
with Eq. (125). 


c Jk+l — (<f>c KccjCk K e I Ik Key l Ik K cx Zk K e $ ft fi c (126) 

Ik+i = h + At(fT y yfc - H z Zk) (127) 

Vxk = K X y yk K xe eh K x j Ik K xz Zk K X f ft — n e (128) 

Of course, it is also necessary to compute and update the command state and forcing 
vectors, Zk and ft, respectively, according to the desired trajectory. Note that since the 
feedforward law only uses the current values Zk and ft, these commands need not be 
available ahead of time and could be real-time pilot inputs or may be computed from 
real-time pilot inputs. 


C. IMPLEMENTATION 

While Eq. (126) - Eq. (128) with the addition of the command state and forcing 
vector constitute the combined feedforward/feedback control law, there sure a number of 
advantages to implementing the control law in incremental form. 

First note that the constant terms n x and /x* have not been determined. Even though 
it is possible to compute these vectors using trim conditions, it is more convenient if the 
control law were not to require these vectors. A second advantage to an incremental 
implementation is the elimination of the integral terms. When the tracking error is, at 
least temporarily, large, the integral state J* can reach very high values. This usually leads 
to the control commands reaching limits. Even though the tracking error may have been 
reduced to small levels, it can take a considerable amount of time before the integral state 
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reaches reasonable levels, and the limiting of the controls is eliminated. The unnecessary 
effects of this phenomenon, referred to as “integral wind-up” can largely be eliminated by 
an incremental implementation. 

By simply differencing the control law, 


AI k = A t(H y y k - 1 - H z z k -\) (129) 

Ack +1 = (<t>c - Kcc) Ac*r- Kcyivie — j/Jfc-i) - Kei AIk 

- Kcz(z k - z k-i) - K es (Ck ~ ?fc-i) (130) 

u x k = - K xy (y k - yk- i) - K xc A c k - K x i A I k 

- K xz (z k - z*-i) — K k ( (ffc - ffc-i) (131) 

z k + 1 = f{ z k >ffc) » Ck = 9k (132) 

When u x k has been modeled as a control rate command with a zero order hold, the actual 
control position commanded is given by 6k. 

6k+ 1 = + At u x k , (133) 

where the control position command 6k is part of the state x*. Other holds will result in 
similar expressions. 

Thus, the actual implemented digital control law is given by Eq. (129) - Eq. (133). 
It is seen that the constant terms depending on the trim conditions have canceled out and 
do not appear in this implementation. 
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It should be noted that the initial condition for the error integrator states, I 0 , is not 
needed in this implementation. Only the initial control variables u I0 , and sometimes S ot 
and the initial compensator increment A c 0 are needed at initialization. The initial condi- 
tions for the control variables can be set equal to the actual control values at initialization. 
When a dynamic compensator is designed, the best selection of the initial compensator 
value is not clear; however, in many cases, the objective of the compensator and the initial 
plant operating point provide a good choice. For example, when the plant initially is in 
trim (i.e., in a steady state condition), the initial compensator increment, Ac 0 , would be 
selected as zero. According to circumstances, other choices are possible. A more detailed 
study of the selection of the initial conditions, particularly when obvious choices are not 
available, is necessary. 

When the plant, due to mechanical reasons, has limiting effects on the movement of 
the actual control variables, it is desirable not to command the controls to exceed these 
limits since such commands will not be followed. Therefore, it is often desirable to have 
control rate and control position limits set in the control law. In this implementation, such 
limits are easy to implement and generally have little negative impact. Rate limits can be 
applied to u x * and position limits to £*. 

It is important to distinguish between limiting action due to the feedforward com- 
mands as opposed to feedback related commands. It is important to note that, in a 
satisfactory design, there should be no plant limiting due to feedforward commands, ex- 
cept in some circumstances. The feedforward control design should include an analysis of 
maximum control commands implied by the class of commanded trajectories. In general, 
the plant limiting conditions should be avoided by appropriate changes in the command 
state and forcing vector. 

In other words, the feedforward design philosophy proposed taken here is to command 
only trajectories which can physically be achieved by the plant, and avoid using up the 
control authority in the feedforward control, thus leaving some control authority to the 
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feedback controller. This has two desirable effects. The first is to allow the feedback 
controller to close the loop and provide its main objective; i.e., stability. This is of utmost 
importance in plants which are open-loop unstable or have relaxed static stability, as is 
the case with many high performance aircraft. If the feedforward commands were to reach 
the plant control limits, the feedback law would not be able to close the loop and perform 
its critical objectives. The second effect of this philosophy is to maintain the nonlinearities 
in the command model generating the desired trajectory. If an unachieveable trajectory is 
commanded, the precise outcome is not clear; i.e., the actual and desired trajectories will 
diverge; however, the nature of the divergence is no longer controlled, and how to recover 
from the divergence is not clear. Whereas by commanding and tracking a trajectory which 
may be Bomewhat different than originally desired, tracking control is maintained, and 
can be used to converge with the originally desired trajectory. The ease with which such 
nonlinear command models can be implemented digitally, as opposed to analog designs, is 
also worthy of note. 

Finally, the usual type of “integrator wind-up” is eliminated in this implementation 
since the integral itself is not explicitly computed. Of course, when no limiting occurs, 
the effect of the integrator is unchanged; the integration is simply performed at a different 
location, namely in Eq. (131). However, when (nonlinear) limiting occurs, the effects are 
usually much more benign. 

Eigenvalues of Implementation. 

Since the implementation is obtained by differencing the control law designed, it would 
result in the same numerical control commands when all initial conditions are appropriately 
matched, no nonlinearities and no random disturbances are present. However, since these 
conditions rarely, if at all, hold, the implemented and designed control commands are not 
the same. 

It is important to note that the implemented control law is closely related to but 
different than the designed law. For example, the implemented law depends on both y* 
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and yjfe-i, whereas the designed law only uses y k . It is therefore necessary to investigate 
the closed-loop characteristics of the implemented control law. 

After some manipulation, the implemented closed-loop system can be written as 

( %k+i \ (<i>x-V x KxyC x — T X K XC —T x K x i T x K xy C x r x \ / x k y 
Ac* + i K cy C x <f> c Kcc K e i K C yC x 0 Ac k 

AI k+1 = A tH y C x 0 0 0 0 A I if (134) 

Xfc I 0 0 0 0 xjf-i 

u xk J V -K xy C x -K xc —Kxi K X yC x iJ \uxk-J 

where the command model state and forcing vector, z k and f*, and the constant vector d k 
have been set to zero as they do not affect or modify the closed-loop eigenvalues. 

It is clear that the stability of the implemented closed-loop system is determined by 
the eigenvalues of the matrix, $ j, rather than the eigenvalues of the designed closed-loop 
system shown below. 

fx k+1 \ (<f>x - TxKxyCx -T X K XC -TxK x \ ( x k \ 

c k+1 = -KcyCx tc-Kcc -K cI C* , (135) 

\l k+1 J \ AtHyCx o I )\l k ) 

where the vectors z k , $ k and the trim related constant vectors, which do not influence the 
closed-loop eigenvalues have again been set to zero for convenience. 

The state transition matrix in Eq. (135) will be denoted by while the implemented 
state transition matrix in Eq. (134) will be denoted by $/. Now, to investigate the 
relationship between the eigenvalues of $/ and $£>, consider 

$/X = AX , X T = (x T c T I T x*u T ) (136) 

where X has n + n s 4- n r + n 4- r components. 

THEOREM 1. Any eigenvalue of $£> is an eigenvalue of $/. 

PROOF: Suppose that ( x T c T J T ) is an eigenvector of corresponding to the eigen- 
value A. Note that there is at least one eigenvector for each eigenvalue, no matter its 
multiplicity. Let 
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x = A® , c = (A — l)c , 7 = (A — 1)7 , x — x , (137) 

u = — K xy C x x- K xc c — K x j 7 , (138) 

for the augmented vector X according to the partition shown in Eq. (136). 

To show that X is eigenvector of corresponding to the eigenvalue A, we write X 
according to the partition and manipulate using the fact that A is also an eigenvalue of 
$D’ 


{4>x - T x K xy C x ){\x) - T X K XC {\ - l)e - T x K xI {\ - 1)7 + T x K xy C x x + T x u 


= A (4> x - T x K xk C x )x - Ar*7f xc c - XT x K xI I 


= A(Ax) = Ax 


(139) 


-K cy C x {Xx) + (< i> c - K CC ){X - l)c - K cI {X - 1)7 + K cy C x x 
= (A - l) (Ac) = A® 


(140) 


AtH y C x {Xx) = A(A7 — 7) = A(A - 1)7 = A7 (141) 


7® = A® = A® 


(142) 


— 7T iy C x (A®) - K xc {X - l)c - K xI (X - 1)7 + K xy C x x + u 
= (A — l)u + u = Au 


(143) 
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Therefore, X is eigenvector of $/ corresponding to the eigenvalue A, and the assertion 
is proved. 

Thus, the designed eigenvalues are maintained in this implementation. Note that we 
have further shown that the eigenvectors of the implementation are very closely related to 
those designed; i.e., Eq. (137), Eq. (138). Note that other incremental implementations 
also maintain the designed eigenvalues [10]. 

It should be noted that the case of eigenvalues with multiplicity greater than one has 
not been considered, although, it would seem that the multiplicity may also be preserved. 

THEOREM 2. has exactly n + ni + n e — r c zero eigenvalues, where 


r c — rank ( <f> e ~ K ce K c j K e yC x ) (144) 

PROOF: Using Eq. (136) for A of zero, we obtain 


<f> x x + r*u = Ax = 0 (145) 

— K cy c x x + ( <f> c — K cc )e — K c jl + K C yC x x _ = Ac = 0 (146) 

A tH y C x x = XI = 0 (147) 

lx = Ax_ = 0 (148) 

u = -K xy C x x - K xc C - K xI I + K xy C x x- + u = Au = 0 (149) 


From Eq. (148) and Eq. (149), it is seen that x and u must vanish. Then, Eq. (145), 
Eq. (147), Eq. (148) and Eq. (149) are automatically satisfied. The remaining constraints 
are 
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{<t> C - K cc )c - Ktfl + K cy C -x_ =0 


(150) 


u = K xc c + K x il — K xy C x x— (X51) 

X r = (0 r C T I T x T _ u T ) (152) 

where u is given by Eq. (151) and (c, I, x_) satisfies Eq. (150) is an eigenvector of 
$/ corresponding to the zero eigenvalue. Since Eq. (150) has n c + ni + n — R c linearly 
independent solutions, there are as many linearly independent vectors X satisfying Eq. 
(150) - Eq. (152), which is the desired result. 

Note that since r c < n c , the implementation has at least n + nj zero eigenvalues. 
When, as is usual, the number of integrals is selected to be the number of plant controls, 
all eigenvalues are accounted by the two theorems. 

COROLLARY. Let $d have distinct non-zero eigenvalues and n j = r; then the eigenvalues 
of $ i consist of those of $£> and zero. 

PROOF: By Theorem 1, the n+n c + ni eigenvalues of $£> also belong to $/. By Theorem 
2, has n + r zero eigenvalues, which are necessarily different than those of $2?* Thus, 
all 2(n + r) + n c eigenvalues of $£> are accounted for. 

It is interesting to note that if the number of integrators used is greater than the 
number of plant controls, then some of the eigenvalues of must necessarily be zero. 
This seems to provide a further implication that the number of integrators should be 
selected to be no greater than the number of plant controls, irrespective of the order of 
the dynamic compensator. 

On the other hand, if the number of integrators is smaller than the number of plant 
controls, all the eigenvalues of are not necessarily accounted for by Theorems 1 and 2. 
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In such cases, it is necessary to investigate the implemented closed-loop eigenvalues further, 
to ascertain that unstable eigenvalues (e.g., A = 1) axe not introduced in the process. To 
accommodate some cases, it is sometimes possible to use different implementations which 
circumvent implementation instabilities [14], 

To illustrate the possibility of unstable eigenvalues, the following Theorem will be 
stated without proof. 

THEOREM 3. $/ has an eigenvalue of unity if, and only if, the matrix below is singular. 

( 4>z- 1 o r x \ 

AtK cI H y C x I-<f> c + K cc 0 (153) 

\AtK xI H y C x K xc 0 J 

This theorem illustrates that when no integrators are used in the design, so that 

- ’ T 

H y vanishes, the matrix in Eq. (153) becomes singular, and the implementation has 
eigenvalues equal to unity which may, and often do, cause problems. This effect of the use 
of integrators provides a further incentive for their use in the control law. 
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IH. DESIGN OF DIGITAL AUTOMATIC LANDING SYSTEM 


In this section, the design of a digital integrated automatic landing system for NASA’s 
Advanced Transport Operating Systems (ATOPS) Research Aircraft will be described. 
The ATOPS aircraft is a Boeing 737-100 which is used by Langley Research Center as a 
research vehicle equipped with special equipment and flight computers with which auto- 
matic control systems can be implemented in flight tests. 

The automatic landing system described in this report is the successor to the Digital 
Integrated Automatic Landing System (DIALS) designed by the author using full state 
feedback techniques and random disturbauice accommodation results [14], [8]. This system 
was successfully flight tested by NASA Langley Research Center [15]. The automatic 
control system described here uses stochastic output feedback [12], [11] and the stochastic 
feedforward techniques developed in the previous section. 

A. Lateral Control Law Design 

The design of the lateral control system follows the basic approach described in the 
previous section; i.e., the stochastic feedforward and output feedback design methodology. 
This requires the development of the design model, the feedback and the feedforward 
controller designs. 

Lateral Plant Design Model. 

As described in the previous section, the design model is needed for both the feed- 
back amd the feedforward control law designs. This model contains the aircraft’s lateral 
aerodynamic and kinematic model, as well as the control actuator models, any inner loop 
control law and filtering which is already present and is intended to be part of the overall 
controller. The design model also includes basic elements of the control law structure such 
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as the open-loop models of any desired dynamic compensation, error integral feedback and 
control rate command. 

The aerodynamic aircraft model can be obtained by linearizing the general nonlinear 
equations of motion (e.g., [16], [17]) about a flight condition with level wings, constant 
airspeed, constant flight path angle corresponding to the desired glideslope, flaps and gear 
down in the landing configuration. The linearized equations obtained are of the following 
form. 


t/ = ant/ + O12 P + &13* + &n<P + a xl W' v 


+ bnSA + b\2$R + bi^Ssp + d v 

(154) 

P = C21V 1 + 0 2 2P + 023»* + <*21^v 

+ b2\6A + & 22 ^-R + b23SSp + dp 

(155) 

r = aait/ + 032 P ■+• 0331 * + <i3iWl 

+ b^iSA + 632^/2 + bzaSsp + d r 

(156) 

= p + tan 0 O r + d+ 

(157) 

ip = r + d^, 

(158) 


where <p , ip, 9 are the roll, yaw and pitch angles, respectively, p and r the roll rate and the 
yaw rate, respectively, v' is the inertial velocity of the aircraft c.g. along the y body-axis 
normalized by the nominal airspeed; 6 A, 6R and 6s p are the aileron, rudder and spoiler 
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control surface positions, respectively; 0 O is the nominal value of the pitch angle. The 
coefficients o,-y and 6,-y are functions of the aircraft’s stability and control derivatives at 
the flight condition used in the linearization. 

The lateral wind model included in the design model is a simple first order system 
describing the normalized lateral wind velocity, iv' u , driven by white noise. 

W' = -0.1W''+u/ w (159) 

The kinematics in the lateral axis consist of the y-position of the aircraft c.g. relative 
to some fixed axis. Figure 2 shows the lateral geometry and the definition of some of the 
relevant singles and distances. The Earth-fixed set of axes has its origin at the glide path 
intercept point (GPIP) , with the positive x-axis pointing along the runway centerline in the 
direction of the aircraft velocity at landing. The z-axis is along the local vertical positive 
downwards; with the y-axis selected so as to make a right-handed coordinate frame. 

Let V„ be the nominal airspeed of the aircraft and Leb — L be the 3x3 matrix 
representing the transformation from the aircraft body axes to the Earth-fixed coordinate 
frame described above. If L,y is the (t, j) element of the matrix L, then 

y' = L 21 u 1 + 1/22 V 1 + Z»23 W* » (1®®) 

where u ; , v' and w' are the normalized inertial velocity components along x, y and z body 
axes, respectively, and 


L 21 = cos 9 sin ip 

(161) 

L 2 2 = sin <p sin 9 sin ip + cos <p cos ip 

(162) 

L 23 = cos ip sin 6 sin ip — sin ip cos ip 

(163) 
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The general expression in Eq. (160) can be approximated using small angle approxi- 
mations to give 


y' = v' - w' a <p + u' cos 0 o %!) + d y , (164) 

where u’ 0 and w' 0 are the normalized nominal velocity components along the x and y body 
axes, respectively. 

The lateral actuator models used for control law design purposes are given below 


6 A = —12.378 6 A + 


(11)(1.531) 
1 + .0015? 


6A C 


(165) 


6R = -32.67 6 R + 32.67(r x + r 2 ) + 32.67 SR C , (166) 


where 6A and SR represent the actual aileron and rudder surface positions, while SA C 
and 6R C represent the aileron and rudder commands generated by the control law, respec- 
tively, and q is the dynamic pressure. The variables ri and r 2 are inner loop (yaw damper) 
variables which will be discussed in more detail later. The aileron and rudder actuator 
models given above are linearized and approximated versions of nonlinear systems contain- 
ing servomechanisms, hydraulic and mechanical systems with usual nonlinearities such as 
hysteresis and limiting effects. A more detailed discussion of the actuator systems on the 
ATOPS research vehicle can be found in [18]. The spoiler is not used as an independent 
control surface, but rather as an aid in producing further rolling moment during turns in 


rAAT^ora f iati wtfV n!lr?. A .7! Q’.l rf "TP 
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on the flight computers, which is approximated to a simple linear relation in the design 
model. 


6s p = 1.73 6 A 


(167) 
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The ATOPS research aircraft incorporates a yaw damper. In the design of the auto- 
matic landing system, the yaw damper is taken to be part of the basic airplane stability 
control system and interpreted as an inner loop system which will be part of the overall 
controller. The yaw damper model given below is therefore included in the design model 
of the open-loop plant. 


fi = -0.3(ri + r 2 ) , 

(168) 

r 2 = -6.993 r 2 + 10.7485 r , 

(169) 

where r is the yaw rate, modeled by Eq. (156). The variables rj and r 2 are then used to 
generate the overall rudder command as shown in the rudder actuator model in Eq. (165). 

Another system that is included in the design model as part of the open-loop plant is 
a third order complementary filter. This filter uses a body-mounted accelerometer triad 
along with position information from the Microwave Landing System (MLS) to obtain 
estimates of the aircraft velocity in the Earth coordinate system. The complementary 
filter is approximated as follows. 

CLy = iio = — lOzio + 10dy , 

(170) 

y ; = i u = Xi 2 + .8(y - Xu) , 

(171) 

y' = x l2 = a'y + .32 (y - i u ) , 

(172) 

a' = ^ = Q + r , 

v u 0 p 

(173) 
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where a v is the body-mounted accelerometer reading for the x-axis, while a' y is its normal- 
ized form and X\q the filtered acceleration, in and x 12 are the filtered aircraft distance 
from the runway centerline and x 12 its rate of change. 

As mentioned earlier, the control variables are the aileron and rudder commands 6A C 
and 6R C , respectively. A rate command structure is used in the design of the control law, 
even though the actual implementation will command the surface positions. 


^•dcfc+i = 6A c k + At uik , 

(174) 

6Rck + 1 = &Rck + At tt2k > 

(175) 


where At is the sampling interval of the control law. 

The design model accommodates a second order dynamic compensator, even though 
the actual design does not use the dynamic compensator. However, the model is given 
here for completeness. 


Cifc+i = « 5At cifc + At u ci k 


(176) 

C2*+l = e -5At C2 k + At u c2 k 

j 

(177) 


Finally, the control law structure is modeled with two integrators, shown below as 
digital accumulators. 


hk+i = Ilk + At(xnk + 5 Xi2k ~ y'ck ~ 5 v'ck ) > 

hk+l — Ilk + A t{<pie — <Pek ) > 

where y c *, y e k and <p c k are the commanded values for y* and <Pk, respectively. These 
commanded values will be defined in the command model. 
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The open-loop plant model used for the design of the lateral control law is obtained 
by the augmented system consisting of the state equations (154) - (159), (164) - (166), 
(168) - (172) and (174) - (179). It is important to note that Eqs. (174) - (179) which 
contain the digital controller structure are defined in discrete form in exactly the way 
that they will be implemented. Whereas the remaining equations of the design model are 
given in differential equation form as they model continuous processes. The latter set of 
equations must first be discretized using the usual sampled-data formulation [19] based on 
the assumption that 6A C and 6R C remain constant over the sampling interval. Then this 
set of discrete equations are augmented by the already discrete set (174) - (179) to obtain 
the complete discrete design model for the open-loop plant of 20 tfc order. 

*fc+i = <t>x%k + r* u* + T, Zfc -|- Wfc + d x (180) 

Lateral Command Design Model. 

While the design model of the open-loop plant developed in the preceding is suffi- 
cient to design the feedback control law, the feedforward control law design requires the 
command design model. The command design model is used to obtain the structure and 
gains of the control law; however, the implementation of the actual commands may use a 
somewhat different set of equations as will be discussed in more detail in the following. 

The lateral command model is selected to be the 4 th order system given by 

A * 2 

y'ck+i = y'ck + At v'ck + -5-y'ck > (isi) 

iick+i = y'ck + A *y'ck , ( 182 ) 


. At 2 . 

Pck + 1 = <Pck + A * <Pck H — <Pck 


(183) 
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£>c*+I = Pck + &t<Pck 


(184) 


The command model state vector, z k , is 


z k = 


[y'ck ii'ck <Pck 



(185) 


Clk 


At 2 


y'ck 


?2 k = At y ck , 


At 2 _ _ 

~Y' (pc k > f4Jfe = Atv?cfc 


The resulting command model is 


(186) 


Zk+l = 4>zZ k + i k 


(187) 


/I 

0 

0 

Vo 


At 0 
1 0 
0 1 
0 0 



(188) 


As can be seen from the open-loop plant design model integral feedback equations 
(178) and (179), the tracking variables are the lateral position y and the roll angle <p. In 
the case of the y-position, a linear combination of the position and velocity is used as 
the error integral feedback. Furthermore, since the actual position and velocity are not 
known, their estimates as obtained from the third order complementary filter are used in 
the feedback loop; thus, the position tracking error is defined to be 


(y'k - y'ck) + 5(yfc - y'ck) > (189) 

while the bank angle error is simply the difference tp k — <p ck . 

An error formulation of the type described in the previous section is used for the 
design of the controller. The variables in which error terms are formed are the roll angle 
<p, the roll rate p, the lateral position y, the yaw angle ip, and the lateral position and 
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velocity estimates from the complementary filter. The yaw angle command used is the 
desired track angle. With this definition of H in the error formulation Eq. (180) becomes 

Xk = Xk- H z k (190) 

S'fc+i = x k 4" r* ttfc 4" (r* + <f > x H H <f > z ) z k H $ k 4* + d x (191) 

The error formulation of the open-loop plant model thus obtained is then used in 
designing the feedforward and feedback control laws. Table 1 shows the feedback vector 
and summarizes the structure of the controller. Table 2 shows the feedforward matrices 
T, and T f which apply to the error formulation of the open-loop design model. 

While the command design model is given by Eqs. (181) - (188), the on-line generation 
of the actually commanded path is obtained using a more complex procedure. The lateral 
trajectory followed by the aircraft may be divided into two portions: the localizer capture 
path and the localizer. The localizer beam is assumed to be on the runway centerline. 

Figure 2 shows the basic geometry of lateral maneuvers. When the controller is 
engaged, the aircraft heading xp a , is extended until it intersects the runway centerline 
(hence the localizer) at xj. It is assumed that the initial aircraft position and heading 
are such that the aircraft would intersect the runway centerline if the heading remained 
constant. A new independent variable, R, is defined as follows. 

R = (x - x/) cos^- + y sin ^ , (192) 

I It 

As can be seen from Figure 2, R is measured along the bisector of the heading t/> 0 at the 
localizer intercept point, x/; it is the position component of the aircraft along this bisector; 
i.e., R is the distance between the localizer intercept point, xj, and the intersection of the 
perpendicular to the bisector. The localizer capture command path is defined with R as 
the independent variable. 
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The command path prior to the initiation of the localizer capture mode simply follows 
a straight line along the initial heading of rp Q . The localizer capture path is commanded 
when the capture criterion is satisfied, at R — R 0 . The localizer capture path is a smooth 
curve with continuous first and second derivatives at both capture initiation and termi- 
nation. At R = R 0 + P, the commanded localizer capture path smoothly transitions into 
the straight-in localizer portion. The decrab mode is initiated when the decrab altitude is 
reached and continues until touchdown. 

Thus, the actual lateral command path is generated using the equations 

(x - x/) tan ip Q , R < R 0 or JLOC = F 

Vc = < fo{R) , Ro < R < Ro + P 

,0 , R > Ro + P 

r Vq sin ip 0 , R < R 0 or JLOC = F 

i/c = < R fi (P) j Ro ^ R < Ro + P 

0 , R > Ro + P 

R = x cos^f- + y sin^- , V G = \Zi 2 + y 2 , (195) 

It Z 

where V G is the estimated ground speed. 

The smoothness of the localizer capture path is due to the selection of the functions 
fo(R) and fi(R). Over the interval (— y , y], these functions are defined by 

fo(R) = fo{Ro) + fi{Ro){R ~ Ro) + Ro) 2 ~ ^ + co 3 ~p~^ ( 196 ) 

fi(R) — fi{Ro) + A (R — Ro) + — sin ^ J , (197) 

h(R)=A(l + cos^ S j , (198) 
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(193) 

(194) 


where fo(Ro), /i(-Ro), Ro, P and A are arbitrary parameters. These parameters are 
selected so as to satisfy the boundary conditions of continuity of the first and second 
derivatives of fo(R), namely fi(R) and / 2 (f2), at both ends of the interval [— y , y] with 
the adjoining paths, resulting in 

p = ~j sin Y . R ° = j ain T * ( 19Q ) 

fo{Ro) = sin 2 y. , / 1 (H 0 )=2«'ny , (200) 

^ , ( 201 ) 

where Vo is the commanded airspeed and |y| max is the maximum inertial lateral accelera- 
tion which will be required to track y e perfectly. Since |y| mo * I s a measure of the sharpness 
of the capture and of the maximum bank angle required, it is left as a parameter to be 
selected by the flight experimenter. When |y| m a* is increased the localizer capture ma- 
neuver will be engaged closer to the localizer and will be performed more quickly using a 
higher bank angle. 

The localizer capture criterion resulting from this trajectory is to engage the capture 
mode when 


ISI ^( 2 V °‘ in T ) 2 < 202 > 

is satisfied for the first time. 

The commands for the roll angle, <p e , are chosen so as to produce a coordinated turn 
when perfect tracking is achieved. The commanded track angle can be found to be 

/i (R) cos 

if> c = tan~ 1 t- , Ro<R<Ro + P . (203) 

1 ~ h{R) sin ~£- 
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Accordingly, the commanded roll is selected to be 



When the decrab altitude is reached, the commanded roll angle is selected so as to 
roll into the wind and perform a sideslip maneuver 

<Pck-G d tp k , <Pck+i = <Pck + . (205) 

The actual command vector used on-line is thus obtained using the set of equations 
described above. On the other hand, the design of the control gains uses the command 
design model given by Equations (181) - (184). Using the design models for the open-loop 
plaint and controller structure and the command trajectory, the feedback and feedforward 
design for the lateral control law is obtained using the stochastic feedforward and output 
feedback approach described in the previous section. The analysis and simulation results 
will be described in the next section. 

B. Longitudinal Control Law Design 

The longitudinal control law design is performed following the same approach as the 
lateral control law. Although the flight maneuvers performed in the longitudinal vertical 
plane sure different than the lateral maneuvers, command path models similar to the lateral 
capture can be used in the glideslope capture and flare maneuvers. As in the lateral case, 
a longitudinal plant design model and a command design model are needed to obtain the 
longitudinal control law. 

Longitudinal Plant Design Model. 

The plant design model is obtained by combining the aircraft’s longitudinal aerody- 
namics and kinematics model with the control actuator and complementary filter models, 
and then augmenting the resulting model with the controller model consisting of dynamic 
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compensation, integral feedback and control rate models. An inner loop controller is not 
included in the longitudinal open-loop plant model. However, dynamic models for the 
vertical position and velocity estimates of the third order complementary filter and for the 
elevator and engine dynamics are included. 

The aircraft’s longitudinal aerodynamic model is obtained by linearizing the nonlinear 
aircraft equations of motion ([16], [17]) about the desired flight condition. The resulting 
equations are of the form 

V' = on V' + a \2 w' + <*13 9 + <*i4 8 + <*ii Wy + <*12 W' w 

+ 6n 8c + 612 87 + dy , (206) 


w' — O21 V* "k <*22 W 1 + 023 9 "k O24 8 + O21 Wy "k <*22 


*k &21 8c ■+■ 622 87 + d w 


(207) 


q = 031 V 1 -k 032 u/ + 033 q + 034 0 + 031 Wy + 032 W w 

• + &31 8 c + 632 87 + dq , (208) 

0 — q + d$ , (209) 

where V 1 is the normalized inertial speed of the aircraft e.g., w' the normalized velocity 
along the body z-axis, q the pitch rate, 0 the pitch angle, W{r and W' w the normalized wind 
velocities along the x and z axes, respectively, 8c the elevator surface position, 8T the engine 
thrust, dy, d w , d p and d e constants depending on the linearization point. The coefficients 
o,-y and bij are constants depending on the aircraft stability and control derivatives [16], 
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[17], [14]. The normalization factor for the affected variables is the aircraft’s nominal 
airspeed, V 0 . 

It should be noted that, due to the presence of the linearization constants dy, d w , d q 
and dg , the state variables V ' , W', q and 0 are not the usual perturbations but the total 
value of these variables. 

The design model for the longitudinal wind velocity components Wy and W' w is taken 
to consist of first order dynamics in each of the velocity components. Thus, 

= -0.1 W{r+w V , (210) 

W , w =- 0 . 1 W , w +u w , (211) 

where «/ and u> w are assumed to be independent white noise processes driving the wind 
model. 

The position of the aircraft c.g. along the Earth-fixed x and z axes can be obtained 
from the kinematic equations. The general kinematic equations can be expressed as 

x' = Vq cos ip T = V' cos 7 cos tpr , (212) 

z' — —V' sin 7 , (213) 

where Vq is the normalized ground speed, tp? the track angle and 7 the flight path angle. 
In the plant design model used here, the position along the Earth-fixed x-axis is taken 
as the independent variable, and is not modeled as a state variable. Accordingly, the 
expression for x given in Eq. (212) is not a state equation. On the other hand, the vertical 
position z is modeled as a state variable. After some manipulation, Eq. (213) can be 
approximated in the form 
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z' = cos 7 „ w' — cos loO — sin 7 „ V' + d z 


(214) 


where is the nominal flight path angle. 

The longitudinal controls used are the elevator surface position, 6e, and the engine 
thrust, ST. The control actuator models used for design purposes are given by 

i n 7 a 

Si = -23.23 Sc + 2.0779 — — — Se c , (215) 

1 + 0.0023? c v ' 

ST = -0.5 ST + 0.298 Sth , (216) 

6th = -Sth + 6th c , (217) 

where Se c and 6th c represent the commanded elevator and throttle positions, respectively, 
and Sth the throttle position. The elevator actuator model includes the effect of cable 
stretch due to aerodynamic loading on the surface. It is also important to note that the 
engine dynamics are actually rather nonlinear, and respond faster when reducing thrust 
than when increasing thrust. A linear approximation to the latter condition has been used 
in the design model in Eq. (216). A 1 second time constant is used to model the throttle 
servo dynamics. Rate limiting and other nonlinear effects are not included in the open-loop 
plant design model. 

The third order complementary filter is modeled by 


z = 0.8s- 0.82 + z 


(218) 


z = 0.4 z — 0.4 z + a z 


(219) 


a z — -10 a z + 10 a z 


( 220 ) 
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a t = -10a x + 10a x ( 221 ) 

where a x and a 2 are the filtered versions of the actual accelerations a x and a z , respectively. 

As in the case of the lateral design model, a second order dynamic compensator is 
included in the plant model; however, the design does not use dynamic compensator and 
purposefully results in vanishingly small gains. The open-loop compensator model is shown 
here for completeness. 


c lfc+l — u clk > (222) 

C2fc+l = ti c 2 k > (223) 

where u e ik and u C 2 * are the compensator control variables. 

To obtain a Type-1 system in the commanded variables, integral feedback of the error 
in altitude and airspeed is included in the control structure as follows. 

hk+i = hk + At(5* + 5 z k — z' ck - 5 z' ck ) , (224) 

/a*+i = hk + At(VZ + WU - V' k ) , (225) 

where At is the sampling interval of the control law. Note that although the inertial 
speed V' is not known, the sum [V' + W' V ) is the airspeed which is measured and therefore 
available for feedback. Vlk here represents the normalized airspeed rather than inertial 
speed command. 

Finally, the design model uses a rate command structure for the control variables. 
Thus, 


Se ck +i = Se ck + At uu 


(226) 
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6th ck+ i = 6th ck + At u 2 * , (227) 

where Ui* and u 2 * are the control variables of the open-loop plant design model thus 
obtained. 

When the continuous system model described by linear differential equations are dis- 
cretized using the standard stochastic sampled-data formulation and the digital controller 
model is added, the discrete plant design model of 20 th order is obtained. 

£k+i = <f>x + T x Ufc + T* Z k + w* + d x . (228) 

To avoid confusing the vertical position variable, z, with the command state vector, 
the latter is denoted by Z in Eq. (228) and in the following. 

Longitudinal Command Design Model. 

The longitudinal variables used in the design, as indicated by the error integral feed- 
back variables in Eqs. (224) and (225), are the airspeed and a linear combination of the 
vertical position and its rate of change. Note that this linear combination may be inter- 
preted as the predicted value of the vertical position in 5 seconds. Also recall that since 
the Earth-fixed z-axis is defined positive downward, the vertical position, z, is the negative 
of the c.g. altitude, h. 

The longitudinal command model selected for designing the feedforward control gains 
is the 2 nd order system given by 


Z\k+x - Z\k + 


(229) 


Z2k+1 = %2 k + k 


(230) 


Zk+i = I Zk + $k 


(231) 
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The command state vector Zk corresponds to ( z ck V c'k) T > where z' ck is the normalized 
vertical position command and V' k is the normalized airspeed command. 

An error formulation is also used in the longitudinal control law. Thus, error terms 
are formed in the normalized vertical position, z, the complementary filtered estimate of 
the vertical position, z, and the airspeed, V. With this definition of the matrix If, the 
plant equations become 


Xk — Xk — H Zk (232) 

S *+1 = <f>x Xk + r x u* + (r* + 4>x H — H<j) z )Zk — H $k + W* + d x . (233) 

While the command design model given by Eqs. (229) and (230) are used to design the 
feedforward control gains, the actually commanded path is obtained as follows. Initially, 
the aircraft is assumed to be in level flight prior to the glideslope capture maneuver. The 
initial altitude of the aircraft is maintained until the glideslope capture criterion is satisfied, 
at which time the desired glideslope capture path is commanded. When the capture is 
completed, the desired glideslope is the vertical path commanded until the flare initiation 
criterion is satisfied. At that time, the altitude profile corresponding to the flare path is 
commanded until touchdown. 

The same functions /o, /i and ji which have been used in the lateral command path 
generation are used in generating the glideslope capture and flare path. An important 
characteristic of these paths is their smoothness. The independent variable used in the 
longitudinal path command is the x-position of the aircraft c.g. in the Earth-fixed coor- 
dinate frame. The commanded path is given here as the altitude, h, which corresponds 
to the negative of the vertical position, z. Thus, the glideslope capture and flare vertical 
profiles are generated using the basic form 
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(234) 


/io(^) = J&o(xo) + /&i(xo)(x — xo) + — A(x — xo) 2 

_ A (^) 2 ( 1+C0S Mf±M) 

hi ( x ) = hi (x 0 ) + A J(x - x 0 ) + ^ sin 27r ( g + Ag ) j f (235) 

hz(x) = .A^l + cos > (236) 

where xo corresponds to the initiation of the maneuver, and the constants are selected 
according to whether the glideslope capture or flare maneuver is to be performed. The 
altitude, sink rate and vertical acceleration are given by 


h(x) = ho(x) , h(x) = xhi(x ) , h(x) = x 2 hi(x) + xhi{x) . (237) 

The glideslope capture path starts from level flight and smoothly transitions to the 
desired glideslope angle ~fcs • When the capture path starts at an altitude hoc and 
requires a vertical acceleration no larger than |h| max , the parameters of the path given by 
Eqs. (234) - (235) are 


Aac = 


|h|ma« 

~2Vf 


Pgc = 


tan igs 

Aqc 


(238) 


hoc Pgc 

Xo = XGC = 7 JT - 

tan igs 2 


Ax = 


hGC 
tan 7 gs 


(239) 


ho(xGc) = hcc » hi(iGc)=0 (240) 

Thus, when the aircraft is at an initial altitude h G c and must track the glideslope 7 G s 
without exceeding |h| ma * of vertical acceleration the vertical profile given by Eq. (234) 
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using the parameters in Eqs. (238) - (240) will be denoted by hac(x ) = h 0 Gc( x ), hiGc{x) 
and hzaci. 1 )- It should be noted that the smoothness of the glideslope capture can be 
easily adjusted by selecting the appropriate value for \h\ max . 

The criterion for initiating the glideslope capture can be expressed as follows. Com- 
mand the path hoGc( x ) when the inequality 


XktaniGS - hk < — tankas (241) 

A 

is satisfied. Here x*, and hk are the current estimates of the aircraft’s x-position and 
altitude. 

For the flare maneuver, the constraints axe placed at touchdown. It is desired that 
the aircraft touch down at xtd with a flight path angle 7 td- Since flare must also initiate 
on the glideslope, the parameters for the flare profile hp(x) are uniquely determined. 


4 _ (tan 'Ud -tan 'las) 2 D _ tan'Ud - tankas 

F 2(ht d — xtdtan'iGs) ’ F A F ' [ } 


Xq = Xp = Xtd — Pf 


A P F 

Ax = — xp 

2 


(243) 


ho(xp) = xptan'iGS > hi (xp) = tankas - (244) 

The flare initiation criterion resulting from this trajectory is to command the flare 
path hp(x) when 


hk 5: xp tan '/as + 12.75 (245) 

is satisfied for the first time. Note that the constant 12.75 allows for the fact that the 
altitude of interest for touchdown is measured to the bottom of the wheels on the landing 
*The glideslope capture criterion was later changed, replacing hk by h„(xk ) to eliminate 
a spike in the elevator. 
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gear rather than the aircraft c.g.. The overall vertical profile of the commanded path 
starting at level flight prior to glideslope capture until touchdown can be expressed as 


h{x) = 


' h-GC 

^Gc(*) 

xtari'ics 


l ^f(z) 


, X < XGC 

i X GC < x < X GS = X GC + Pgc 
, X GS < x < X F 
, X > Xp 


( 246 ) 


M*) = { 


h 2 {x) = 


(o 

hlGc[x) 
tan igs 

l hipfe) 

( 0 

h 2 Gc{ x ) 

0 

l h 2 f (*) 


, X < x G c 

, xgc < x < xgs = xgc + Pgc 
, %gs < x < x f 
, x > x f 

i 

, x < xgc 

, x gc < x < x gs = x gc + Pgc 

, X GS < x < 

, X > X/p 


( 247 ) 


( 248 ) 


Recall that the sink rate and vertical acceleration can be easily obtained using Eq. 


( 237 ). 

The airspeed command is generated as follows. Let V Q be the commanded airspeed, 
V a the initial measured airspeed and V c k the airspeed commanded at the k th sample by the 
command model. Since a sudden jump in airspeed is undesirable, any difference between 
the initial and desired airspeeds, V Q - V 0 , is gradually eliminated starting with an initial 
command equal to the actually measured airspeed. 


Vdk+1 = V ck + fcfc 


V' = V' 

Y CO r o 


V ck = 


Vck 


( 249 ) 


?2 k = < 


sgn(V 0 - 


. if |l — ^cfcl > 


( 250 ) 


^ 1 — V* k , otherwise 

where V max is the maximum acceleration or deceleration desired during the approach. 
Thus, the commanded airspeed starts from the actual one and commands a constant 
deceleration or acceleration until the desired airspeed command is reached. 
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During flare, the airspeed command decelerates the aircraft for touchdown as long as 
the airspeed remains above the minimum desirable airspeed V m i n which is selected higher 
than the stall speed and consistent with contingencies such as go-around maneuvers. 


v; k+l 


-{ 


K* + ft* 


nun 


Vik > VL 


ck - 


vi, 


(251) 


ftfc+1 = S 


& k - At 


25 


V 0 Tn 


-At 


25 

V a T 


lo 


, f 2 * +1 > -At — and VSk > 

, tok+i < -At ^ and VS k > Cin 

, VSk < C 


(252) 


where T is the period of time in which a 25 ft/sec. decrease in airspeed will be commanded. 

Thus, during flare the commanded airspeed is reduced to aid in the maneuver and in 
touchdown as long as a safe speed of V m i n is maintained. 
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IV. ANALYSIS AND NONLINEAR SIMULATION 


The lateral and longitudinal/vertical control laws designed use the open- loop plant 
design models and the command design models described in the last section. The design 
approach used is the Stochastic Feedforward/Output Feedback methodology described in 
Section II. In this section, the design is analyzed by obtaining the closed-loop eigenvalues 
and singular values, and by simulating the digital automatic landing system obtained to 
control a nonlinear computer simulation of the ATOPS Research Vehicle, a B-737-100 
aircraft. 


A. CLOSED- LOOP SYSTEM ANALYSIS 

The open-loop plant design models for the lateral and longitudinal/vertical dynamics 
have been discussed in Section III. As can be seen by simple observation of the equations 
making up the design models, both the lateral and longitudinal/ vertical open-loop plants 
can be expressed in the form 

Xfe+i = <j>x xic + r* u* + r , Zk + r f & + w* + d x (253) 

The lateral design model state consists of the twenty components: 


{v' t P, r, <P, y', W ' v , S A, SR, a' y , y\ r x , r 2 , c u c 2 , I\, / 2 , 6 A c , 6R e } 

The control law structure is selected so that 13 out of the 20 lateral states are used 
in the feedback. In particular, the actual aileron and rudder surface positions, the yaw 
damper inner loop system states and the lateral wind velocity are not used for feedback, 
although some of these are measured and available for use. The feedback vector for lateral 
controller consists of 
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{P,r,<p,ii>,a , v ,y',y',ci,c 2 ,l u Ii,6A e ,6R e } 

The control vector consists of the aileron and rudder rate commands SA C and SR C , 
respectively. Integral feedback is used for errors in lateral position and roll angle. The lat- 
eral control structure is summarized in Table 1. Table 2 summarizes the lateral command 
model parameters needed in designing the feedforward control gains. Table 3 shows the 
designed lateral gains for the feedforward and feedback control laws. The control gains in 
Table 3 follow the terminology of Eq. (128). 

It should be noted that the feedforward command state gain corresponds to the case 
of error feedback. This is obtained intentionally to aid in steady state offset reduction. The 
reasoning can be illustrated as follows. Suppose that the altitude and sink rate errors and 
the airspeed error are null,, and that the command system does not contemplate a maneu- 
ver, then it is reasonable to maintain the control surfaces at their current position, since 
changing the control commands to different values is likely to result in a non-zero error. 
The same approach is used in both the lateral and longitudinal feedforward controllers. 

Closing the loop with the feedback control gains given in Table 3 results in the closed- 
loop discrete system. The s-domain equivalent eigenvalues of the closed-loop lateral system 
are shown in Table 4. The singular values for the discrete plant with the loop broken at 
the input is shown in Table 5. 

The longitudinal/vertical design model state consists of the following twenty compo- 
nents: 


{V\w,q,9y t Wl t Wl,6T t 6e,8thXX'*\2*u°*'Ii*l2,6th et 6e e } 


The longitudinal feedback vector excludes some of the components of the state vector 
such as the inertial speed in vertical body axis u> # , the longitudinal and vertical wind 
velocities, the true aircraft altitude, the engine thrust and the elevator surface position. 


68 


On the other hand, the complementary filtered altitude (—5), sink rate and accelerations 
are used by the controller, as can be seen by the feedback vector components 

{V r/ ,9,tf,a , x ,o , y ,2 / ,i',c i,c 2 ,h,l2,6th e ,6e c } 

In the longitudinal/ vertical control law, the commanded variables are predicted vertical 
position and airspeed. Thus, the integral of the error in the predicted vertical position 
and the airspeed are used in the feedback law. The throttle rate and elevator rate are, the 
control components; however, as discussed in Section II C, the actual control commands 
are the throttle amd elevator positions 6th c and 6e Ct respectively. The longitudinal/vertical 
feedback control structure is summarized in Table 6. 

The command design model parameters required for the feedforward control of the 
longitudinal/vertical trajectory is summarized in Table 7. The gains designed for the 
feedforward and feedback controllers are shown in Table 8. 

Closing the loop with the output feedback gains obtained results in the closed-loop 
equivalent s-domain eigenvalues shown in Table 9. The singular value, eigenvalue and Bode 
plots of the closed-loop system sure shown in Figure 4. 

B. NONLINEAR SIMULATION 

The performance of the digital automatic landing system described above is evaluated 
in this part through a digital computer simulation. The ATOPS B-737-100 aerodynamics, 
actuator systems, kinematics, servo, hydraulic and other systems have been simulated in 
considerable detail in a nonlinear digital computer simulation. In this simulation, dynamic 
systems such as the complementary filter, the yaw damper, the spoiler-aileron coupling, 
the engine, etc. are modeled as nonlinear systems which accurately describe their actual 
behavior rather than their linearized versions used in the open-loop plant model. 

The digital automatic control system described in the preceding sections is simulated 
in detail. The control law simulation is then interfaced with the aircraft simulation so that 
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the control commands computed by the design become the input to the aircraft control 
actuator systems. Numerous simulations of the closed-loop aircraft system were performed 
under a variety of different conditions. The aircraft response in simulations are shown in 
Figures 5 - 16. 

The digital control system simulates the incremental implementation shown in Section 
II. C Eqs. (129) - (133). While other incremental implementations are possible, they 
are not used here. It should be noted that, since both lateral and longitudinal/ vertical 
controllers designed used a control rate structure, Eq. (133) is, in fact, implemented to 
obtain the commands for the aileron, rudder, elevator and throttle positions. The actual 
outputs of the control system are the control position commands 6 *. 

The Boeing 737 aircraft used here has a baseline stabilizer automatic trim logic. 
The auto-trim logic drives the stabilizer surface so as to minimize the moment on the 
elevator hinge, thus providing maximum authority for the elevator to react to sudden 
changes in the flight parameters. The stabilizer movement is much slower than that of 
the elevator and does not introduce further dynamic modes in the models. The use of the 
incremental implementation is very suitable to accommodate such slow moving surfaces. 
In the nonlinear simulations of the automatic landing system designed here, the stabilizer 
automatic trim logic is turned on; so that the stabilizer automatically trims the aircraft 
even though the plo.ts shown do not include the stabilizer position. 

The control system iteration rate used in the simulation is 10 Hz which is also the 
sampling rate used in the control law design. Since the control law is digital, this update 
rate must be used for the controller simulation. However, the simulation of the aircraft 
aerodynamic and on-board systems is performed at 20 Hz. Since these system describe 
continuous processes with some modes of high natural frequency, their simulation requires 
a higher update rate for accuracy. It is also important to note from Eqs. (129) - (133) 
that the control system output £* at time f* uses only variables available at time j. 
Therefore, as long as the real-time computation of the commands, require no more 
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than 100 msecs, no computational delay will be present. It is assumed that a sufficiently 
fast flight computer will be used to compute the incremental implementation commands. 
Thus, the accuracy of the digital control system simulation is expected to be restricted 
to round-off errors due to limitations in the word length of the flight computer, possible 
mismatches in obtaining exactly a 10 Hz rate, input-output limitations, etc. 

The feedforward and feedback control gains used in the simulation are given in Tables 
3 and 8. The simulation of the actual commanded path shown by Eq. (132) in the 
incremental implementation equations is performed as described in Section III. It should 
be noted that the actual command model uses estimates of the position and of the velocity 
of the aircraft c.g.. thus, the feedforward control law implementation actually contains 
nonlinear feedback. The coupling is usually rather low and may be neglected. However, a 
more complete evaluation should include these feedback effects as well as the effect of the 
feedback control law. 

The simulations shown in Figures 5-16 are initialized at an estimated altitude of 
950 ft. At initialization, the aircraft is flying a constant altitude path with level wings 
and a heading so as to intercept the localizer or the runway centerline at some point, xj. 
The automatic landing system is engaged at the initialization of the simulation. At that 
point, the control law checks to see whether the localizer glideslope capture criteria are 
satisfied. The initial conditions have been selected so that neither the localizer nor the 
glideslope capture criterion will be met at this point. Thus, according to the command 
path generated, the aircraft continues along the same track angle with level wings and 
maintains a constant altitude. 

At initialization, the aircraft calibrated airspeed is selected to be 135 knots. On the 
other hand, in most of the simulations, the commanded airspeed V a is 125 knots. Accord- 
ingly, the control law experiences an instantaneous error of about 10 knots at initialization. 
As described in Section III, the feedforward command model generates a linearly decreasing 
airspeed profile from this initial speed to its commanded speed. As seen in the simulations, 
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the aircraft decelerates, following the commanded airspeed, until V a is reached. When the 
commanded airspeed V 0 is 135 knots, as shown in Figure 6, then the initial airspeed error 
is zero, and the airspeed command is a constant value of V 0 . As can bee seen from Figure 
6, the control law accurately maintains the commanded speed of 135 knots rather than 
reducing it to 125 knots as it does when the commanded airspeed is 125 knots. 

As the control law continually tests to see if the glideslope or localizer capture criteria 
are satisfied, whether the glideslope or the localizer capture will be the first to engage 
depends on the aircraft’s position, its heading relative to the runway centerline and the 
commanded glideslope angle, with other parameters having generally a lesser effect. In 
almost all the simulations shown here, one capture mode is engaged soon after the other, 
so that the aircraft flight path is a curved 3D path in both the lateral and vertical planes. 
Ability to perform both localizer and glideslope captures simultaneously is described in 
order to achieve close-in captures, as it is no longer necessary to perform localizer capture 
first and then engage the glideslope capture. 

When the desired glideslope angle is 3°, the glideslope capture criterion is satisfied 
first in most of the simulations. The initiation of the glideslope capture maneuvers can be 
clearly seen in the commanded and actual sink rate plots. Both commanded and actual 
sink rate smoothly transition from level flight to the sink rate required to remain on the 
desired glideslope at the desired airspeed. Also note the pitch angle and angle of attack 
movements during the glideslope capture maneuvers. Whereas the angles coincide when 
flying a constant altitude path, when glideslope capture is engaged automatically, the 
control law smoothly pitches the aircraft down to capture the glideslope. Note that there 
is no initial tendency to pitch in the “wrong” direction. 

Also note that prior to glideslope capture, the pitch and angle of attack have to 
increase slightly when the aircraft is decelerating in order to compensate for the reduction 
of lift due to the airspeed, hence dynamic pressure, reduction. The needed extra lift is 
obtained by pitching up and increasing the angle of attack, albeit with lag which results 
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in a small altitude offset. On the other hand, when the aircraft does not decelerate, as 
shown in Figure 6, the aircraft pitch and angle of attack remain essentially constant as it 
is not necessary to obtain extra lift to maintain altitude. 

It should be noted that since maneuvers are principally performed by the feedfor- 
ward controller, the glideslope capture performance is an indication that the feedforward 
controller is satisfactory for this type of maneuver. Observation of the throttle and eleva- 
tor positions shows that the feedforward controller pitches down by initially lowering the 
throttle rather than using the elevator. For the B-737 which has a thrust line considerably 
lower than the c.g., reducing thrust has the added effect of reducing the pitching moment. 
This control strategy is precisely the one that best suits this aircraft, since simply using 
the elevator to pitch would have the unwanted result of increasing the airspeed. Thus, the 
stochastic feedforward control approach is indeed making use of the plant design model 
information as would be desired. 

At the end of the glideslope capture, the altitude error is redefined causing an instan- 
taneous move in its value. This causes a corresponding sudden and undesired transient 
in the elevator and throttle commands when the altitude error move is appreciable. This 
glitch can be removed in a number of ways, including the use of a simple easy-on function 
at the appropriate time. 

It should be noted that the smoothness of the glideslope capture maneuver and its du- 
ration are directly related to the parameter |/i| mai ; i.e., the maximum vertical acceleration 
of the commanded altitude profile. By simply varying this parameter at any time prior to 
glideslope capture, the commanded capture path may be changed on-line to a smoother 
or faster maneuver as desired. 

In Figures 7 and 8, the automatic landing system captures steep glideslope of 4° and 
4.5°, respectively. The glideslope capture path generated by the feedforward model given 
in the previous section by Eqs. (234) - (248) is automatically modified to result in a vertical 
profile which captures a steeper glideslope, tracks it and flares from this steeper glideslope. 
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Since the same maximum vertical acceleration value is used for all glideslope angles, it may 
be noted that the duration of the steeper glideslope captures is slightly longer. However, 
the basic characteristics of the response and of the control law remain unchanged. 

It should also be noted that the feedforward and feedback control gains remain un- 
changed while the commanded vertical path varies according to speed or glideslope angle. 
Thus, the feedforward control law tracks the commanded path within the class of com- 
manded trajectories. It is clear that if the commanded paths are sufficiently different 
from each other, feedforward controllers adapted to the specific characteristics of each 
path would result in “better” performance. One such approach would be to extend the 
stochastic feedforward approach to include optimal gain scheduling. 

In the simulation of the steep glideslope cases, it may be noted that glideslope capture 
occurs later. This is a consequence of the geometry, as the initial aircraft altitude is the 
same. In the case of the 4° glideslope both localizer and glideslope capture occur in the 
same period of time. Localizer capture occurs mostly before the aircraft captures the 
steeper 4.5° glideslope in part due to a larger localizer intercept angle of 47°. In all cases, 
the control law captures the desired glideslope by satisfactorily tracking the commanded 
capture path. 

When the glideslope capture is performed, the aircraft tracks the desired glideslope 
until the flare mode js engaged. As can be seen from the simulations of different glideslope 
angles (i.e., Figures 5-8), the aircraft remains on the desired glideslope with essentially 
the same precision for shallow and steep glideslopes. 

The flare mode is engaged when the flare criterion in Eq. (245) is satisfied. In this 
mode, both the flare vertical profile generated as well as the airspeed reduction profile are 
commanded. As can be seen from the simulations, the aircraft pitches up increasing its 
angle of attack and the lift as desired; this results in a corresponding reduction in the sink 
rate and the airspeed until touchdown. This is achieved by using the elevator to pitch up 
while lowering the throttle to reduce the airspeed. 
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It can be seen that in all cases, the pitch angle at touchdown is comfortable above 
zero and still rising. This pitch attitude at touchdown is necessary to avoid landing on the 
nose wheel which is not designed for the high load at touchdown. 

As the flare maneuver is significantly sharper than the glideslope capture maneuver, 
it requires faster control action. Thus, the measurement noise covariance, V" f , for the 
command model forcing vector was reduced to obtain the flare gains, resulting in a higher 
forcing vector gain, K s , during flare as shown in Table 8. Due to the complexity of the 
flare maneuver and the high accuracy needed in tracking the altitude profile, a higher order 
altitude command model would model this profile more accurately. 

The localizer capture for a 3° glideslope and localizer intercept angle of 32°, as shown 
in Figure 5, is initiated near the end of the glideslope capture maneuver. As can be seen 
from the heading and roll plots, the aircraft yaw and roll angles track their commanded 
trajectories closely and capture the localizer. The lateral position is also seen to track its 
commanded profile accurately. Although a small deviation from the localizer is present, 
this does not resemble a usual overshoot pattern as it occurs after reaching the localizer. 
In all cases, this offset is quite small and tends not to exceed 3 m. As in the case of 
the glideslope capture, the high accuracy of the tracking of the lateral position indicates 
that the stochastic feedforward controller can produce satisfactory feedforward control law 
designs. 

During the maneuver, the roll angle and the lateral position commands are selected 
so as to produce a coordinated turn when perfect tracking occurs. The sideslip angle plots 
indicate that the sideslip angle remains well within 1° of sideslip during the whole final 
approach excluding, of course, the decrab maneuver shown in Figure 14. The maximum 
sideslip tends to occur during localizer capture slightly after the peak bank angle. On the 
other hand, note that the lateral acceleration is plotted on the same set of axes and follows 
the roll angle quite closely, as expected during a coordinated turn, whereas the lateral 
specific force in the body axes remains near zero in a coordinated turn. 
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The smoothness and duration of the localizer capture path command can be varied 
by varying the maximum commanded lateral acceleration |y| mai . This maximum acceler- 
ation is set to 10 ft/sec 2 in most of the simulations shown. However, Figure 10 shows a 
simulation where |y| ma x is set to a value of 5 ft/ sec 2 . Observation of the lateral variables 
shows that the localizer capture is smoother, lasts longer, requires a lower maximum roll 
angle to capture as well as resulting in a lower maximum lateral acceleration as expected. 
The turn coordination, although acceptable for |y|ma* value of 10, appears to be improved 
as indicated by a lower maximum sideslip angle. The tracking of the lateral position also 
shows some improvement. Thus, the feedforward control law can clearly track lateral com- 
mand paths where the lateral acceleration does not exceed 10 ft/ sec 2 with no change in 
the feedforward control gains. 

A number of simulations are shown in Figures 5-16 where the sensitivity of the 
automatic landing system to various parameters is shown. For example, the aircraft weight 
is varied among 85,000 lbs, 90,000 lbs, and 95,000 lbs. The center of gravity of the aircraft 
is also varied in tandem with the weight. Figures 13 - 16 shows the sensitivity to winds 
and noise including bias errors. The wind gust standard deviation used in the simulations 
containing gusts is 2 ft/sec. (0.61 m/sec.). The airspeed command is varied between 125 
and 135 knots, while the commanded glideslope angles simulated are 3°, 4° and 4.5°. 
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V. CONCLUSIONS AND RECOMMENDATIONS 


In this study, a combined stochastic feedforward /feedback control design methodology 
is developed, and a digital automatic landing system is designed using this approach. It 
is considered that the main objective of a control law is to enable the plant to track a 
desired or commanded trajectory selected from a given class of trajectories as closely as 
possible in the presence of random and deterministic disturbances and despite uncertainties 
about the plant. The feedforward controller tries to track the desired or commanded 
trajectory, whereas the feedback controller tries to maintain the plant state near the desired 
trajectory despite the presence of random, and possibly deterministic, disturbances and 
uncertainties about the plant. Modern control theory has concentrated more attention 
on the important feedback control problem, while the feedforward control problem has 
received less attention. 

The feedforward control problem is formulated as a stochastic output feedback prob- 
lem where the plant contains unstable and uncontrollable modes. As the standard output 
feedback algorithm requires an initial gain which stabilizes the plant, a new algorithm is 
developed to obtain the feedforward control gains. The necessary conditions are shown to 
result in coupled linear matrix equations, implying that when a solution exists, it is indeed 
the globally optimal control gain. 

The formulation of the feedforward problem in a stochastic, rather than the standard 
deterministic, setting is significant in two ways. First, the class of desired trajectories 
from which the actually commanded path is selected can be effectively described as a 
random process generated by a dynamical system driven by a white noise process. The 
second, and more important, implication of a stochastic optimization formulation is the 
tacit understanding that “perfect tracking” is often not possible due to various reasons 
including uncertainties about, or variation in the, plant parameters, the presence of plant 
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nonlinearities and unmatched initial conditions. Thus, questions about the robustness and 
sensitivity of the feedforward controller arise naturally in this context. 

A combined stochastic feedforward/feedback control methodology is developed where 
the main objectives of the feedforward and feedback control laws are clearly seen. Fur- 
thermore, the inclusion of error integral feedback, dynamic compensation, rate command 
control structure, etc. is an integral element of the methodology. Another advantage 
of the methodology is the flexibility that a variety of feedback control design techniques 
with arbitrary structures may be employed to obtain the feedback controller; these include 
stochastic output feedback, multi-configuration control, decentralized control or frequency 
and classical control methods. 

Finally, a specific incremental implementation is recommended for the combined feed- 
forward/feedback controller. Some advantages of this digital implementation are the sim- 
plicity of implementation, the fact that trim values are not needed and that problems 
such as integrator wind-up can be largely avoided. The closed-loop eigenvalues using this 
implementation are shown to contain the designed closed-loop eigenvalues which would re- 
sult if an incremental implementation were not used. It is further shown that when using 
an incremental implementation, it is advantageous to design the controller with as many 
integrators as the number of controls. Using fewer integrators results in marginally stable 
eigenvalues of unity, while using more integrators constrains the placement of eigenvalues. 
The choice of the same number of integrators as controls is also an intuitively pleasing one. 

A digital automatic landing system for the ATOPS Research Vehicle (a Boeing 737- 
100) is designed using the stochastic feedforward controller and stochastic output feedback. 
The system control modes include localizer and glideslope capture, localizer and glideslope 
track, crab, decrab and flare. Using the recommended incremental implementation, the 
control laws are simulated on a digital computer and interfaced with a nonlinear digital 
simulation of the aircraft and its systems. 

In this study, the feedforward controller takes an equal place along the feedback con- 
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troller in achieving the overall control objective. While the stochastic feedforward/feedback 
approach has been successfully developed and applied to a significant problem, some signif- 
icant questions and extensions of the problem remain unanswered, and are recommended 
for further study and experimentation. Three general areas of study are worthy of further 
investigation: 

• the structure of the feedforward controller 

• the robustness and sensitivity of the feedforward controller 

• optimal gain scheduling of the feedforward controller 

The structure of the feedback controller considers questions about the role of feedfor- 
ward dynamic compensation, the use of the “future values of the desired trajectory” in the 
current control command, the use of the full-state feedforward controller when fast flight 
computers are available. An argument can be effectively made that since a pilot knows 
the future desired trajectory and uses this information in his current control commands, 
the optimal feedforward controller should also take advantage of such information. 

The uncertainties about complex system parameters and nonlinear effects bring forth 
uncertainties about the trajectory which would be tracked when the actual plant parameter 
are different than those used in the feedforward design. Since the feedforward controller 
does not determine the stability of the closed-loop plant, instability does not generally 
result from such mismatching. However, since unsatisfactory performance would generally 
result from a high sensitive feedforward law, it is of interest to study measures of robustness 
and design methods which incorporate low sensitivity criteria. 

In applications where the plant will vary over a wide range of conditions resulting in 
large changes in plant model parameters, or in cases where the command model parameter 
vary to achieve some objective, it is necessary to adapt the feedforward control gains 
according to varying conditions. This can be achieved by extending the optimal gain 
scheduling studies to include feedforward controller. Due to the relative simplicity of the 
coupled linear necessary conditions, gain scheduling with respect to all the plant parameters 
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rather than a selected few may be feasible. In particular, the feedforward gain of the 
command model forcing vector seems extremely appropriate for such application. 
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